Childhood exposure to non-persistent endocrine disruptors, glucocorticosteroids, and attentional function: A study based on the parametric g-formula

Abstract

Evidence suggests that endocrine disrupting chemicals (EDCs) may perturb the hypothalamic-pituitary-adrenocortical (HPA) axis, which has a major role in brain development. We aimed to evaluate the effects of childhood exposure to organophosphate pesticides, phenols, and phthalate metabolites, on urinary glucocorticosteroids and inattention in children using data from the Human Early-Life Exposome (HELIX) cohort. We used the parametric g-formula to estimate effects between EDCs, glucocorticosteroids, and hit reaction time standard error (HRT-SE), a measure of inattention from the Attention Network Test (ANT), and tested for possible effect modification by sex. We observed a positive marginal contrast (MC) for exposure increases from the 10th to the 90th percentile for methyl-paraben (MC: 0.042 and \(95\%\) confidence interval (CI): (0.013, 0.071)), and the phthalate metabolites oxo-MiNP (MC: 0.023 and \(95\%\) CI: (0.003, 0.044)), oh-MiNP (MC: 0.039 and \(95\%\) CI: (0.001, 0.076)), and MEHP (MC: 0.036 and \(95\%\) CI: (0.008, 0.063)), on HRT-SE, indicating lower attention. Several EDCs were also associated with a positive MC for cortisone, cortisol, and corticosterone production. Increased levels of the glucocorticosteroids had no effect on HRT-SE, although we found a possible effect modification by sex. Our results suggest that multiple EDCs might interfere with inattention in children and with the homeostasis of the HPA axis.

The prevalence of several neurodevelopmental disorders has increased in the pediatric population (1), and multiple environmental pollutants may play a role in the increased rates of these disorders (2). Multiple endocrine disrupting chemicals (EDCs), ubiquitous chemicals present in many every-day products and diet, are capable of interfering with the endocrine system, and have shown associations with childhood neurodevelopment and behavior (317). Although both pregnancy and early infancy are crucial stages of (neuro)development, most of the available literature is focused on the effects of prenatal exposure to EDCs on child neurodevelopment (2).

One group of EDCs that may have a deleterious effect on neurodevelopment is the organophosphate pesticides (OP pesticides), although the few studies assessing exposure during childhood and through the use of biomarkers suffered from a series of limitations, including a small sample size (2). Exposure to phthalates and their metabolites during childhood and early adolescence has also been associated with several adverse neurodevelopmental outcomes, but these studies were limited to few phthalate metabolites and small study populations (2). The effects of exposure to bisphenol A (BPA) during childhood on cognitive functions are still unclear (2).

Moreover, little is known about the biological mechanisms of action (2). There is some toxicological evidence, however, that exposure to certain EDCs, specifically phthalates, might interfere with the hypothalamic-pituitary-adrenocortical (HPA) axis and might interact with the glucocorticoid receptor (1820). The HPA axis, which can be activated by stress, is responsible for the production of glucocorticosteroids. The brain, and its proper functioning, is a potential target, due to the presence of receptors for these hormones (19,21). Glucocorticosteroids are necessary for brain maturation, although their under- or over-production might interfere with its normal development and ultimately lead to long-term impaired functioning (20,21).

Taken together, these results suggest that the negative influence of exposure to certain EDCs on neurodevelopmental outcomes might be mediated, at least partially, by disruption of the HPA axis’ homeostasis. In the present study, we thus estimated cross-sectional associations between 1) non-persistent EDCs and attentional function, 2) non-persistent EDCs and glucocorticosteroids, and 3) glucocorticosteroids and attentional function, using the parametric g-formula and marginal contrasts (MCs), in children of a large network of cohorts in Europe.

Methods

Study population and design

The Human Early-Life Exposome (HELIX) project aims to characterize early-life exposures and their potential association with endogenous biomarkers and health outcomes (22). It consists of six existing population-based birth cohort studies across Europe: BiB (Born in Bradford, UK) (23), EDEN (Study of determinants of pre- and postnatal developmental, France) (24), INMA (Environment and Childhood, Spain) (25), KANC (Kaunas Cohort, Lithuania) (26), MoBa (The Norwegian Mother and Child Cohort Study, Norway) (27), and Rhea (Mother–Child Cohort in Crete, Greece) (28). The HELIX subcohort of 1,301 mother-child pairs was fully characterized for the external and internal exposome, including exposure and omics biomarkers during childhood (29). Eligibility criteria for inclusion in the HELIX subcohort included: a) age 6-11 years, with a preference for 7-9 years; b) availability of sufficient stored pregnancy blood and urine samples; c) availability of complete address history from first to last follow-up; d) no serious health problems, which might affect the results of the clinical testing. Ethical permission was obtained from the relevant authorities in the corresponding country.

Variables

Endocrine disrupting chemicals

Children were assessed between December 2013 and February 2016, and assessments included neurological testing and urine collection. Urine samples of the night before and the first morning void on the day of the visit were combined to provide a more reliable exposure assessment. Non-persistent EDCs assessed in the urine samples included phthalate metabolites, phenols, and organophosphate (OP) pesticide metabolites. A list of the environmental chemicals determined in urine samples and used for the present study is given in Table S1. Briefly, we analyzed a total of 7 phenols (bisphenol A (BPA), ethyl-paraben (ETPA), methyl-paraben (MEPA), n‑butyl‑paraben (BUPA), oxybenzone (OXBE), propyl-paraben (PRPA), triclosan (TRCS)), 6 non-specific organophosphate pesticide metabolites (diethyl dithiophosphate (DEDTP), diethyl phosphate (DEP), diethyl thiophosphate (DETP), dimethyl dithiophosphate (DMDTP), dimethyl phosphate (DMP), dimethyl thiophosphate (DMTP)), and 10 phthalate metabolites (mono benzyl phthalate (MBzP), monoethyl phthalate (MEP), mono‑2‑ethyl 5‑carboxypentyl phthalate (MECPP), mono‑2‑ethylhexyl phthalate (MEHP), mono‑2‑ethyl‑5‑hydroxyhexyl phthalate (MEHHP), mono‑2‑ethyl‑5‑oxohexyl phthalate (MEOHP), mono‑4‑methyl‑7‑hydroxyoctyl phthalate (oh-MiNP), mono‑4‑methyl‑7‑oxooctyl phthalate (oxo-MiNP), mono‑iso‑butyl phthalate (MiBP), mono‑n‑butyl phthalate (MnBP)) originating from 6 distinct phthalate parent compounds. The laboratory protocols for the analysis are described elsewhere (30).

Glucocorticosteroids

Urine samples of the night before the day of the visit were used to measure levels of the glucocorticosteroids. These included glucocorticosteroids, glucocorticosteroid metabolites, glucocorticosteroid precursors, glucocorticosteroid precursor metabolites, androgens, and androgen metabolites. A list of the glucocorticosteroids determined in urine samples and used for the present study is given in Table S2.

To assess the levels of glucocorticosteroids and their metabolites, LC-MS/MS analysis was applied at the Applied Metabolomics Research Group, IMIM (Hospital del Mar Medical Research Institute). The laboratory protocols for the analysis are described elsewhere (31,32).

Three additional markers, total cortisol production, total cortisone production, and total corticosterone production, were computed based on the following: cortisol production as the sum of cortisol and its metabolites (20α-dihydrocortisol (20aDHF), 20β-dihydrocortisol (20bDHF), 5α,20α-cortol (5a20acortol), 5α,20β-cortol (5a20bcortol), 5α-tetrahydrocortisol (5aTHF), 5β,20α-cortol (5b20acortol), 5β,20β-cortol (5b20bcortol), 5β-dihydrocortisol (5bDHF), 5β-tetrahydrocortisol (5bTHF), 6β-hydroxycortisol (6OHF)), cortisone production as the sum of cortisone and its metabolites (20α-dihydrocortisone (20aDHE), 20β-dihydrocortisone (20bDHE), 5α-tetrahydrocortisone (5aTHE), 5β,20α-cortolone (5b20acortolone), 5β,20β-cortolone (5b20bcortolone), 5β-tetrahydrocortisone (5bTHE), 6β-hydroxycortisone (6OHE)), and corticosterone production as the sum of 11-dehydrocorticosterone (A), 17-deoxycortolone (17-DO-cortolone), 5α-tetrahydrocorticosterone (5aTHB), 5β-tetrahydrocorticosterone (5bTHB).

Attentional function

Cognitive and motor function outcomes were assessed with standardized, non-linguistic, and culturally blind computer tests, including the Attention Network Test (ANT) (33), which provides a measure of efficiency of attentional function. The tests were administered in a standardized way, and with minimal interference from the field workers. Further information can be found in (29). The outcome of interest for the present study is the hit reaction time standard error (HRT-SE) (34), a measure of response speed consistency throughout the test. A high HRT-SE indicates highly variable reaction times, and is considered a measure of inattentiveness.

Confounders

For each research question, defined by a specific type of exposure and outcome, the minimal set of covariates for inclusion in the analyses was selected on the basis of a directed acyclic graph (DAG) built with DAGitty (35) and ggdag (36). The sets of covariates were selected to estimate the total effect of the exposure on the outcome. For effect estimation of the EDCs on glucocorticosteroids and of glucocorticosteroids on HRT-SE, these sets were also sufficient to estimate direct effects. Sample-specific creatinine values were used to adjust for possible dilution effects. Further, each minimal adjustment set was augmented with precision covariates, defined as the set of parents variable of the outcome that are not parents of the exposure. Common confounders were cohort, ethnicity, sex, age, height, weight, and head circumference of the child, consumption of fish, fruit, vegetables, organic food, anf fast food, maternal tobacco consumption, family financial situation and affluence scale (FAS). Models for estimating the effects of EDCs on HRT-SE were further adjusted for child breastfeeding, prenatal maternal active and passive smoking, urine creatinine, child mood and rest before assessment, child neuropsychological diagnosis, marital status, season, and fasting time before assessment. Models for estimating the effects of EDCs on glucocorticosteroids were further adjusted for urine creatinine, season, and fasting time before assessment. Models for estimating the effects of glucocorticosteroids on HRT-SE were further adjusted for child breastfeeding, prenatal maternal active and passive smoking, marital status, EDCs, urine creatinine, child mood and rest before assessment, and child neuropsychological diagnosis. The adjustment sets are provided in the Supplementary Material as text files compatible with DAGitty. Codebooks for the used covariates, by research question, are provided in Supplementary Tables 3, 4, 5.

Statistical methods

Data pre-processing

Concentrations of the glucocorticosteroids were classified as quantifiable, below the limit of quantification (LOQ), possible interference or out of range, and not detected. For each metabolite, we computed the fraction of values below the LOQ and not detected, both within each cohort and overall. We proceeded to impute these values using half the value of the corresponding LOQ, for those metabolites that had less than 30% of missings within each cohort and 20% of missings overall. Information about the lower limit of quantification (LLOQ) for the glucocorticosteroids is provided in Table S6. The remaining missing values were imputed using kNN from the VIM R package (37), for those metabolites that had less than 40% of remaining missings within each cohort and 30% of remaining missings overall. We used 5 nearest neighbors. We natural log-transformed them to improve model fit, assessed with posterior predictive checks. To do so, replicated data were simulated with the fitted models and compared to the observed data. We used the check_predictions function from the performance R package using the default arguments (38). Values of total cortisol, cortisone, and corticosterone production were expressed in nanograms per millilitre (ng/ml).

Concentrations of the non-persistent EDCs were classified as quantifiable, below the limit of detection (LOD), possible interference or out of range, and not analysed. Concentrations below the LOD were singly imputed using a quantile regression approach for the imputation of left-censored missing data, as implemented in the impute.QRILC function from the imputeLCMD R package (39). Information about the lower limits of detection can be found in (30). Chemicals with more than 70% of observations below the LOD were excluded from the present study. Remaining missing values were imputed similarly using kNN. Values of the chemicals were expressed in \(\mu\)grams per litre (\(\mu\)g/L).

Missing values in the clinical outcome were imputed similarly using kNN. We natural log-transformed these to improve model fit, assessed with posterior predictive checks. Values of the clinical outcome were expressed in milliseconds (ms).

Missing values in the covariates were imputed similarly using kNN. Categorical covariates were imputed using the maxCat function, which chooses the level with the most occurrences. Creatinine values were expressed in grams per litre (g/L).

Estimation of balancing weights

To reduce the effect of measured confounders on the exposure-outcome association, stabilized balancing weights were estimated using the energy method available in the WeightIt R package (40). This method estimates weights by minimizing an energy statistic related to covariate balance (41), thus avoiding the need to specify a parametric model. Weights below the 0.1 and above the 0.9 quantiles were trimmed. Trimming might lead to decreased covariate balance and potentially change the estimand, but can also decrease the variability of the weights. Covariate balance was assessed using functionalities provided by the cobalt R package (42). Specifically, we used Love plots to visualize covariate balance before and after adjusting.

G-computation

We estimated MCs with the parametric g-formula, a method of standardization. The parametric g-formula involves the following steps: 1) fit a outcome model including both covariates and balancing weights; 2) create two new datasets identical to the original one but with the exposure shifted according to a user-specified intervention set by a deterministic function of the observed exposure levels; 3) use the outcome model to compute adjusted predictions in the two counterfactual datasets; 4) compute the difference between the means of the adjusted predictions in the counterfactual datasets. The causal parameter of interest was thus specified as the difference in the expected counterfactual outcomes under the shifted exposure levels \(\left( \mathbb{E} \left[ Y^{d_1} \right] - \mathbb{E} \left[ Y^{d_2} \right] \right)\). In order for this parameter to be identified, the usual causal identifiability conditions (no unmeasured confounding, positivity, and consistency) are required. Since these conditions are likely not satisfied, we focused on the estimation of a statistical estimand that is as close as possible to the causal parameter of interest.

We fit the outcome model using the glm function and a Gaussian family with identity link from base R. The exposure variable was modeled using natural cubic splines with 3 degrees of freedom, to more flexibly capture the average dose-response function (ADRF).

To estimate the MCs, we used the avg_comparisons function from the marginaleffects R package (43). The two counterfactual datasets were obtained by setting the exposures levels to 90th percentile (\(d_1\)) and the 10th percentile (\(d_2\)), for each cohort separately. The MCs were computed using the estimated balancing weights above. Robust standard errors were computed with the sandwich R package, using cohort as variable indicating clustering of observations (44,45). For each outcome, we report the results as differences between MCs.

The R code to reproduce analyses and results is available online (https://github.com/lorenzoFabbri/paper-helixSC-neuro).

Effect-modification analysis

We further estimated separate MCs for possible effect-modification by sex. To do so, balancing weights were estimated separately for each level of the sex variable, and an interaction term between the exposure and sex was included in the outcome model. Similarly, the MCs were aggregated separately for each level of sex.

Results

Table 1 and Table S7 provide descriptive statistics for the outcome and covariates for the HELIX subcohort and for each cohort, respectively. Of the 1,301 children of the HELIX subcohort, 1,297 had measurements of the non-persistent EDCs. Measurements of the glucocorticosteroids were available for 1,004 children, of which 980 were matched to the HELIX subcohort. Measurements of both non-persistent EDCs and glucocorticosteroids were available for 976 children of the subcohort. A flowchart describing the sample size for each research question is presented in Figure S1. The sample consisted of 55% males. The median HRT-SE was 300 ms (interquartile range (IQR), 231-368), with lower median values for EDEN, MOBA, and INMA, corresponding to the cohorts with older children. At the time of visit, the median age of the children was 8.06 years. The children were mostly Caucasian (90%), and the largest minority were of Pakistani origin (6.2%).

Levels of unprocessed non-persistent EDCs, after imputation of values below the LOD, and glucocorticosteroids, are presented in Table 2, Table 3, and Table S8. Supplementary Figures 2 and 3 provide information on the measurement classification of the EDCs and glucocorticosteroids by cohort, respectively.

The effective sample sizes before and after balancing weights estimation are presented in Supplementary Tables 9, 10, 11, while basic summary statistics of the estimated balancing weights are presented in Supplementary Tables 12, 13, 14. As expected, the median value of the weights for each exposure was close to \(1.00\).

Figure 1 presents the forest plot for the MCs on the logarithmic scale of the non-persistent EDCs on HRT-SE. For most EDCs, a cohort-specific increase in the levels of the exposures from the 10th to the 90th percentiles was associated with a positive MC, indicating an increase in the values of HRT-SE and thus lower attention. Most of the confidence intervals (CIs) included the null effect, though. Significant effects were observed for the paraben MEPA (MC: 0.042 and \(95\%\) CI: (0.013, 0.071)), and the phthalate metabolites oxo-MiNP (MC: 0.023 and \(95\%\) CI: (0.003, 0.044)), oh-MiNP (MC: 0.039 and \(95\%\) CI: (0.001, 0.076)), and MEHP (MC: 0.036 and \(95\%\) CI: (0.008, 0.063)). The organophosphate pesticide (OP pesticide) DETP was negatively associated with HRT-SE (MC: -0.026 and \(95\%\) CI: (-0.054, 0.001)).

Figure 2 presents the forest plot for the MCs on the logarithmic scale of the non-persistent EDCs on total cortisone, cortisol, and corticosterone production. For most EDCs, a cohort-specific increase in the levels of the exposures from the 10th to the 90th percentiles was associated with a positive MC, indicating an increase in the total production of these metabolites. Exceptions were BUPA, which was associated with negative MCs for all three outcomes, and MiBP, which was associated with a negative MC for total cortisone production only. The majority of the effects for the phenols and phthalate metabolites included the null. The phenol BPA showed the largest MCs across all three outcomes (cortisone production, MC: 0.263 and \(95\%\) CI: (0.131, 0.394); cortisol production, MC: 0.274 and \(95\%\) CI: (0.107, 0.441); corticosterone production, MC: 0.285 and \(95\%\) CI: (0.106, 0.464)).

Figure 3 presents the forest plot for the MCs on the logarithmic scale of the glucocorticosteroids on HRT-SE. All MCs included the null, with no clear indication of directionality of the effect.

Effect modification by sex

Basic summary statistics of the estimated balancing weights for effect modification are presented in Supplementary Tables 15, 16, 17. As expected, the median value of the weights for each exposure was close to \(1.00\).

Table 4 presents the results of the difference between estimates of the MCs on the logarithmic scale for females and males, for the EDCs on the glucocorticosteroids and HRT-SE. For HRT-SE, significant differences were present for the phenol OXBE (MC: 0.032 and \(95\%\) CI: (0.004, 0.061)) and the phthalate metabolites MEP (MC: -0.053 and \(95\%\) CI: (-0.138, 0.033)) and MbZP (MC: -0.066 and \(95\%\) CI: (-0.126, -0.007)). For the glucocorticosteroids, significant differences were present across all three classes of EDCs and for all outcomes. The largest differences were attributable to the OP pesticides DMTP (cortisol production, MC: -0.21 and \(95\%\) CI: (-0.326, -0.094)) and DETP (corticosterone production, (MC: -0.16 and \(95\%\) CI: (-0.332, 0.011)); cortisone production, (MC: -0.097 and \(95\%\) CI: (-0.269, 0.076))). The forest plots of the individual MCs are presented in Supplementary Figures 4 and 5.

Table 5 presents the results of the difference between estimates of the MCs on the logarithmic scale for females and males, for the glucocorticosteroids on HRT-SE. Significant differences were present for cortisone production (MC: 0.14 and \(95\%\) CI: (0.019, 0.261)) and corticosterone production (MC: 0.126 and \(95\%\) CI: (0.009, 0.243)). Furthermore, for all exposures, the MCs had opposite sign (positive for males and negative for females). The forest plot of the individual MCs is presented in Figure S6.

Discussion

The impact of exposure to EDCs on human health has attracted considerable research interest. While research in this area has mainly investigated the effects of prenatal exposure on child neurodevelopment (2), little is still known about childhood exposure. In this study, consisting of 1,297 children from 6 European birth cohorts, we observed that short-term childhood exposure to certain non-persistent EDCs was associated with attentional function (MEPA, MEHP, oh-MiNP, and oxo-MiNP), and with total production of cortisol, cortisone, and corticosterone (DEP, DMP, DMTP, BPA, ETPA, MEPA, MEHP, oh-MiNP, and oxo-MiNP). Increased production of these glucocorticosteroids did not seem to affect attentional function. Some of these effects differed for females and males, including significant differences for the effects of increased production of cortisone and corticosterone on HRT-SE. Specifically, an increased production of these glucocorticosteroids was associated with lower values of HRT-SE for females, and higher values for males. Taken together, these results suggest that these non-persistent EDCs might be responsible for perturbations of the HPA axis’ homeostasis, and that higher levels of these glucocorticosteroids might interfere with different functions of attention in a sex-specific manner.

To the best of our knowledge, no other study has investigated the effects of childhood exposure to multiple classes of non-persistent EDCs in relation to attentional function. More generally, the literature on non-persistent EDCs and neurodevelopment in children has mostly focused on OP pesticides (3,4,6,8), phthalate metabolites (5,9,10,15,17,4648), and BPA (7,13,14). González-Alzaga et al. and Cartier et al. evaluated cross-sectional associations between dialkylphosphate (DAP) metabolites and subtests of the Wechsler Intelligence Scale for Children (49) in European children with ages between 6 and 11 years. Higher levels of DAP metabolites (DMP, DMTP, DMDTP, DEP, DETP, and DEDTP) were associated with lower scores of intelligence quotient (IQ) and verbal comprehension, especially in boys (4), while higher levels of diethylphosphate metabolites (DEP, DETP, DEDTP) were associated with lower working memory scores (6). There is also preliminary evidence of a possible association between exposure to certain OP pesticides and Attention-Deficit / Hyperactivity Disorder (ADHD) in children (3,8). Specifically, Bouchard et al. found evidence of a cross-sectional association between dimethyl alkylphosphate metabolites (DMP, DMTP, and DMDTP) and ADHD in children aged 8 to 15 years from National Health and Nutrition Examination Survey (NHANES), while Yu et al. found a dose-response relationship between DMP and ADHD in Taiwanese children aged 4 to 15 years. Preliminary evidence is also available for several phthalate metabolites in relation to cognitive development in childhood. Higher levels of di(2-ethylhexyl) phthalate metabolites (including MEHP, MEHHP, and MEOHP) were associated with lower intelligence scores in children aged 2 to 12 years (5), lower scores of IQ and verbal intelligence, more omission errors (a measure of inattention), and higher scores of response time variability (a measure of sustained attention) in 6-year old Korean children (10), poorer fine motor skills in preadolescent boys (47), and lower intelligence scores in 7-year old children (17). Further associations were found for MEOHP with lower scores of IQ (5) and verbal intelligence in Taiwanese children aged 6 to 12 years (9), and for dibutyl phthalate metabolites (MnBP and MiBP) with impaired verbal intelligence (9). There is further preliminary evidence that associations between certain phthalate metabolites and cognitive abilities vary by timing of exposure assessment (46). Among phenols, some studies provide preliminary evidence of an association between BPA and ADHD in children aged 8 to 15 years (7) and in a case-control study of children aged 6 to 12 years (13), especially in boys. Except for working memory, there does not seem to be evidence of an association between BPA and cognitive abilities in Spanish boys aged 9 to 11 years (14). Few studies have looked into different classes of non-persistent EDCs. Shoaff et al., for instance, investigated cross-sectional associations between multiple EDCs and ADHD-related behaviors in 15-year old adolescents, finding a higher risk of ADHD-related behavior problems with higher levels of antiandrogenic phthalate metabolites (MEHP, MEHHP, MEOHP, MECPP, MnBP, MiBP, MBzP, monohydroxyisobutyl phthalate (MHiBP), monocarboxyoctyl phthalate (MCOP), monoisononyl phthalate (MNP), and monohydroxybutyl phthalate (MHBP)), especially in boys (15).

We are not aware of other epidemiological studies investigating childhood exposure to phthalates metabolites, phenols, and OP pesticides, in relation to urinary glucocorticosteroid levels in childhood. Prior epidemiological research provides preliminary evidence for an association between certain non-persistent EDCs with higher levels of glucocorticoids (1820). Repeated measures up to 15 months of age of the phthalate metabolites MEHHP, MEOHP, MiBP, and MnBP showed positive associations with free cortisol in Korean children, with no effect modification by sex (18). In a cohort of Chinese pregnant women, phthalate metabolites were measured at 14, 24, and 36 weeks of gestation, and the glucocorticoids cortisol and cortisone were measured in cord blood. Third-trimester levels of MEHP were positively associated with cortisol, while MECPP and MEOHP were negatively associated with cortisone (19). Time- and chemical-dependent sex differences were also found: during the third trimester, MEHHP and MEOHP were positively associated with cortisol in females, while negatively associated in males (19). In a longitudinal study, a mixture of several phthalate metabolites, driven by MEP, MiBP, and MBzP, measured in childhood, showed a positive association with hair cortisol measured at 12 years of age (20). While in the present study we did find positive MCs between some phthalate metabolites (MEHP, oh-MiNP, and oxo-MiNP) and the glucocorticosteroids, there are important differences with the previous studies. First, exposure assessment was performed during gestation (19) or the first 15 months of life (18), not during childhood. Second, the glucocorticosteroids were measured in other matrices, specifically in cord blood (19) or hair (20). Finally, (20) investigated mixture effects. Contrary to these studies (18,20), we did find effect modification by sex.

Adding to these epidemiological studies, previous toxicological research provide evidence for the inhibition by phthalates of human 11\(\beta\)-hydroxysteroid dehydrogenase 2 (11\(\beta\)-HSD2) activity, responsible for the conversion of active cortisol into inactive cortisone (50,51). There is also in silico evidence suggesting that BPA, a phenol, and Triazophos (TAP), a organophosphorus insecticide, can bind to the human glucocorticoid receptor (52,53).

We are also not aware of prior epidemiological studies specifically investigating the effects of elevated levels of glucocorticosteroids in relation to attentional function, although there is evidence that under- or over-production of glucocorticosteroids interfere with the normal development of the brain (21). While we did find sex-specific evidence of an effect, their clinical relevance is questionable.

Our findings should be interpreted in light of the following limitations and strengths. Limitations include the cross-sectional design of the present study. Importantly, the non-persistent EDCs were measured in a pool of night and morning urine samples before the clinical visit, to represent exposure over the previous day, whereas the glucocorticosteroids were measured in the night urine sample. Although we included a wide range of confounders there is the possibility, as with other observational studies, of residual confounding, which might lead to a bias away from the null. Some of the confounders indicated in the adjustment sets had to be removed due to large fractions of missing values. There is further the possibility of misspecification of the outcome model, although we included a spline of the exposure to relax some of the linearity assumptions. The use of more data-adaptive learners was excluded due to the relatively small sample size. We finally acknowledge the possibility that some of chemicals might not act independently (mixture effect). Further research is thus warranted.

Strengths of the present study include the use of pooled urine samples for chemical assessment to obtain more representative long-term exposures, since it is known that these specific EDCs have very short half-lives (54,55). We decided to model both the treatment mechanisms, for the estimation of balancing weights, and the outcomes, with traditional covariates adjustment, to try to obtain doubly robust effect estimates. Finally, we decided not to interpret our results by focusing on the estimated coefficients of possibly misspecified regression models, but by making use of the g-computation procedure and estimate MCs.

In conclusion, in a study of 1,297 children from 6 European birth cohorts, we observed that (i) exposure to non-persistent EDCs in childhood might have short-term effects on HRT-SE in childhood, (ii) exposure to non-persistent EDCs in childhood might disrupt the HPA axis in childhood, and (iii) disruption of the HPA axis in childhood might have short-term, sex-specific effects on HRT-SE. Future studies should investigate how glucocorticosteroids might mediate the adverse effects of exposure to non-persistent EDCs on childhood neurodevelopment (too broad) in larger populations.

References

1.
Grandjean P, Landrigan PJ. Neurobehavioural effects of developmental toxicity. Lancet Neurol. 2014;13(3):330–338.
2.
Ramírez V, Gálvez-Ontiveros Y, González-Domenech PJ, et al. Role of endocrine disrupting chemicals in children’s neurodevelopment. Environmental Research [electronic article]. 2022;203:111890. (https://www.sciencedirect.com/science/article/pii/S0013935121011853). (Accessed January 31, 2023)
3.
Bouchard MF, Bellinger DC, Wright RO, et al. Attention-Deficit/Hyperactivity Disorder and Urinary Metabolites of Organophosphate Pesticides. Pediatrics [electronic article]. 2010;125(6):e1270–e1277. (https://doi.org/10.1542/peds.2009-3058). (Accessed December 29, 2023)
4.
González-Alzaga B, Hernández AF, Rodríguez-Barranco M, et al. Pre- and postnatal exposures to pesticides and neurodevelopmental effects in children living in agricultural communities from South-Eastern Spain. Environment International [electronic article]. 2015;85:229–237. (https://www.sciencedirect.com/science/article/pii/S0160412015300593). (Accessed December 29, 2023)
5.
Huang H-B, Chen H-Y, Su P-H, et al. Fetal and Childhood Exposure to Phthalate Diesters and Cognitive Function in Children Up to 12 Years of Age: Taiwanese Maternal and Infant Cohort Study. PLOS ONE [electronic article]. 2015;10(6):e0131910. (https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0131910). (Accessed December 29, 2023)
6.
Cartier C, Warembourg C, Le Maner-Idrissi G, et al. Organophosphate Insecticide Metabolites in Prenatal and Childhood Urine Samples and Intelligence Scores at 6 Years of Age: Results from the MotherChild PELAGIE Cohort (France). Environmental Health Perspectives [electronic article]. 2016;124(5):674–680. (https://ehp.niehs.nih.gov/doi/10.1289/ehp.1409472). (Accessed December 29, 2023)
7.
Tewar S, Auinger P, Braun JM, et al. Association of Bisphenol A exposure and Attention-Deficit/Hyperactivity Disorder in a national sample of U.S. children. Environmental Research [electronic article]. 2016;150:112–118. (https://www.sciencedirect.com/science/article/pii/S0013935116302110). (Accessed December 29, 2023)
8.
Yu C-J, Du J-C, Chiou H-C, et al. Increased risk of attention-deficit/hyperactivity disorder associated with exposure to organophosphate pesticide in Taiwanese children. Andrology [electronic article]. 2016;4(4):695–705. (https://onlinelibrary.wiley.com/doi/abs/10.1111/andr.12183). (Accessed December 29, 2023)
9.
Huang P-C, Tsai C-H, Chen C-C, et al. Intellectual evaluation of children exposed to phthalate-tainted products after the 2011 Taiwan phthalate episode. Environmental Research [electronic article]. 2017;156:158–166. (https://www.sciencedirect.com/science/article/pii/S0013935116309112). (Accessed December 29, 2023)
10.
Kim JI, Hong Y-C, Shin CH, et al. The effects of maternal and children phthalate exposure on the neurocognitive function of 6-year-old children. Environmental Research [electronic article]. 2017;156:519–525. (https://www.sciencedirect.com/science/article/pii/S0013935116312828). (Accessed December 29, 2023)
11.
Furlong MA, Herring A, Buckley JP, et al. Prenatal exposure to organophosphorus pesticides and childhood neurodevelopmental phenotypes. Environmental Research [electronic article]. 2017;158:737–747. (https://www.sciencedirect.com/science/article/pii/S0013935117301044). (Accessed January 31, 2023)
12.
Braun JM. Early-life exposure to EDCs: Role in childhood obesity and neurodevelopment. Nat Rev Endocrinol [electronic article]. 2017;13(3, 3):161–173. (https://www.nature.com/articles/nrendo.2016.186). (Accessed December 19, 2023)
13.
Li Y, Zhang H, Kuang H, et al. Relationship between bisphenol A exposure and attention-deficit/ hyperactivity disorder: A case-control study for primary school children in Guangzhou, China. Environmental Pollution [electronic article]. 2018;235:141–149. (https://www.sciencedirect.com/science/article/pii/S0269749117327884). (Accessed December 29, 2023)
14.
Rodríguez-Carrillo A, Mustieles V, Pérez-Lobato R, et al. Bisphenol A and cognitive function in school-age boys: Is BPA predominantly related to behavior? NeuroToxicology [electronic article]. 2019;74:162–171. (https://www.sciencedirect.com/science/article/pii/S0161813X19300543). (Accessed December 29, 2023)
15.
Shoaff JR, Coull B, Weuve J, et al. Association of Exposure to Endocrine-Disrupting Chemicals During Adolescence With Attention-Deficit/Hyperactivity DisorderRelated Behaviors. JAMA Network Open [electronic article]. 2020;3(8):e2015041. (https://doi.org/10.1001/jamanetworkopen.2020.15041). (Accessed December 29, 2023)
16.
17.
Vilmand M, Beck IH, Bilenberg N, et al. Prenatal and current phthalate exposure and cognitive development in 7-year-old children from the Odense child cohort. Neurotoxicology and Teratology [electronic article]. 2023;96:107161. (https://www.sciencedirect.com/science/article/pii/S0892036223000119). (Accessed November 16, 2023)
18.
Kim JH, Lee J, Moon H-B, et al. Association of phthalate exposures with urinary free cortisol and 8-hydroxy-2’-deoxyguanosine in early childhood. Science of The Total Environment [electronic article]. 2018;627:506–513. (https://www.sciencedirect.com/science/article/pii/S0048969718301475). (Accessed December 26, 2023)
19.
Sun X, Li J, Jin S, et al. Associations between repeated measures of maternal urinary phthalate metabolites during pregnancy and cord blood glucocorticoids. Environment International [electronic article]. 2018;121:471–479. (https://www.sciencedirect.com/science/article/pii/S0160412018315113). (Accessed November 9, 2023)
20.
Sears CG, Liu Y, Lanphear BP, et al. Evaluating mixtures of urinary phthalate metabolites and serum per-/polyfluoroalkyl substances in relation to adolescent hair cortisol: The HOME Study. American Journal of Epidemiology [electronic article]. 2023;kwad198. (https://doi.org/10.1093/aje/kwad198). (Accessed November 16, 2023)
21.
Lupien SJ, McEwen BS, Gunnar MR, et al. Effects of stress throughout the lifespan on the brain, behaviour and cognition. Nat Rev Neurosci [electronic article]. 2009;10(6, 6):434–445. (https://www.nature.com/articles/nrn2639). (Accessed December 31, 2023)
22.
Vrijheid M, Slama R, Robinson O, et al. The human early-life exposome (HELIX): Project rationale and design. Environ Health Perspect. 2014;122(6):535–544.
23.
Wright J, Small N, Raynor P, et al. Cohort Profile: The Born in Bradford multi-ethnic family cohort study. International Journal of Epidemiology [electronic article]. 2013;42(4):978–991. (https://doi.org/10.1093/ije/dys112). (Accessed December 19, 2023)
24.
Heude B, Forhan A, Slama R, et al. Cohort Profile: The EDEN mother-child cohort on the prenatal and early postnatal determinants of child health and development. International Journal of Epidemiology [electronic article]. 2016;45(2):353–363. (https://doi.org/10.1093/ije/dyv151). (Accessed November 13, 2023)
25.
Guxens M, Ballester F, Espada M, et al. Cohort Profile: The INMAINfancia y Medio Ambiente—(Environment and Childhood) Project. International Journal of Epidemiology [electronic article]. 2012;41(4):930–940. (https://doi.org/10.1093/ije/dyr054). (Accessed November 13, 2023)
26.
Grazuleviciene R, Danileviciute A, Nadisauskiene R, et al. Maternal Smoking,GSTM1 and GSTT1 Polymorphism and Susceptibility to Adverse Pregnancy Outcomes. International Journal of Environmental Research and Public Health [electronic article]. 2009;6(3, 3):1282–1297. (https://www.mdpi.com/1660-4601/6/3/1282). (Accessed November 13, 2023)
27.
Magnus P, Irgens LM, Haug K, et al. Cohort profile: The Norwegian Mother and Child Cohort Study (MoBa). International Journal of Epidemiology [electronic article]. 2006;35(5):1146–1150. (https://doi.org/10.1093/ije/dyl170). (Accessed November 13, 2023)
28.
Chatzi L, Plana E, Daraki V, et al. Metabolic Syndrome in Early Pregnancy and Risk of Preterm Birth. American Journal of Epidemiology [electronic article]. 2009;170(7):829–836. (https://doi.org/10.1093/aje/kwp211). (Accessed November 13, 2023)
29.
Maitre L, Bont J de, Casas M, et al. Human Early Life Exposome (HELIX) study: A European population-based exposome cohort. BMJ Open [electronic article]. 2018;8(9):e021311. (https://bmjopen.bmj.com/content/8/9/e021311). (Accessed November 13, 2023)
30.
Haug LS, Sakhi AK, Cequier E, et al. In-utero and childhood chemical exposome in six European mother-child cohorts. Environment International [electronic article]. 2018;121:751–763. (https://www.sciencedirect.com/science/article/pii/S016041201831225X). (Accessed January 17, 2022)
31.
Marcos J, Renau N, Casals G, et al. Investigation of endogenous corticosteroids profiles in human urine based on liquid chromatography tandem mass spectrometry. Analytica Chimica Acta [electronic article]. 2014;812:92–104. (https://www.sciencedirect.com/science/article/pii/S0003267013015791). (Accessed November 13, 2023)
32.
Gomez-Gomez A, Pozo OJ. Determination of steroid profile in hair by liquid chromatography tandem mass spectrometry. Journal of Chromatography A [electronic article]. 2020;1624:461179. (https://www.sciencedirect.com/science/article/pii/S0021967320304313). (Accessed November 13, 2023)
33.
Rueda MR, Fan J, McCandliss BD, et al. Development of attentional networks in childhood. Neuropsychologia [electronic article]. 2004;42(8):1029–1040. (https://www.sciencedirect.com/science/article/pii/S0028393204000041). (Accessed April 12, 2023)
34.
Sunyer J, Esnaola M, Alvarez-Pedrerol M, et al. Association between Traffic-Related Air Pollution in Schools and Cognitive Development in Primary School Children: A Prospective Cohort Study. PLOS Medicine [electronic article]. 2015;12(3):e1001792. (https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.1001792). (Accessed November 21, 2023)
35.
Textor J, van der Zander B, Gilthorpe MS, et al. Robust causal inference using directed acyclic graphs: The R package “dagitty.” International Journal of Epidemiology [electronic article]. 2016;45(6):1887–1894. (https://doi.org/10.1093/ije/dyw341). (Accessed September 22, 2023)
36.
Barrett M. Ggdag: Analyze and Create Elegant Directed Acyclic Graphs. 2023;(https://github.com/r-causal/ggdag)
37.
Kowarik A, Templ M. Imputation with the R Package VIM. Journal of Statistical Software [electronic article]. 2016;74:1–16. (https://doi.org/10.18637/jss.v074.i07). (Accessed January 17, 2022)
38.
Lüdecke D, Ben-Shachar MS, Patil I, et al. performance: An R package for assessment, comparison and testing of statistical models. Journal of Open Source Software. 2021;6(60):3139.
39.
Lazar C. imputeLCMD: A collection of methods for left-censored missing data imputation. R package, version. 2015;2.
40.
Greifer N. Cobalt: Covariate balance tables and plots. 2023.
41.
Huling JD, Greifer N, Chen G. Independence Weights for Causal Inference with Continuous Treatments. Journal of the American Statistical Association [electronic article]. 2023;0(0):1–14. (https://doi.org/10.1080/01621459.2023.2213485). (Accessed September 15, 2023)
42.
Greifer N. WeightIt: Weighting for covariate balance in observational studies. 2023.
43.
Arel-Bundock V. Marginaleffects: Predictions, comparisons, slopes, marginal means, and hypothesis tests. 2023.(https://marginaleffects.com/)
44.
Zeileis A. Econometric computing with HC and HAC covariance matrix estimators. Journal of Statistical Software. 2004;11(10):1–17.
45.
Zeileis A, Köll S, Graham N. Various versatile variances: An object-oriented implementation of clustered covariances in R. Journal of Statistical Software. 2020;95(1):1–36.
46.
Li N, Papandonatos GD, Calafat AM, et al. Identifying periods of susceptibility to the impact of phthalates on children’s cognitive abilities. Environmental Research [electronic article]. 2019;172:604–614. (https://www.sciencedirect.com/science/article/pii/S0013935119301379). (Accessed March 13, 2024)
47.
Balalian AA, Whyatt RM, Liu X, et al. Prenatal and childhood exposure to phthalates and motor skills at age 11 years. Environmental Research [electronic article]. 2019;171:416–427. (https://www.sciencedirect.com/science/article/pii/S0013935119300507). (Accessed March 13, 2024)
48.
Jankowska A, Polańska K, Hanke W, et al. Prenatal and early postnatal phthalate exposure and child neurodevelopment at age of 7 years – Polish Mother and Child Cohort. Environmental Research [electronic article]. 2019;177:108626. (https://www.sciencedirect.com/science/article/pii/S0013935119304232). (Accessed March 13, 2024)
49.
Wechsler D, Kodama H. Wechsler intelligence scale for children. Psychological corporation New York; 2003.
50.
Zhao B, Chu Y, Huang Y, et al. Structure-dependent inhibition of human and rat 11β-hydroxysteroid dehydrogenase 2 activities by phthalates. Chemico-Biological Interactions [electronic article]. 2010;183(1):79–84. (https://www.sciencedirect.com/science/article/pii/S0009279709003950). (Accessed February 13, 2024)
51.
Ma X, Lian Q-Q, Dong Q, et al. Environmental inhibitors of 11β-hydroxysteroid dehydrogenase type 2. Toxicology [electronic article]. 2011;285(3):83–89. (https://www.sciencedirect.com/science/article/pii/S0300483X11001466). (Accessed February 13, 2024)
52.
Prasanth GK, Divya LM, Sadasivan C. Bisphenol-A can bind to human glucocorticoid receptor as an agonist: An in silico study. Journal of Applied Toxicology [electronic article]. 2010;30(8):769–774. (https://onlinelibrary.wiley.com/doi/abs/10.1002/jat.1570). (Accessed February 13, 2024)
53.
Yang F-W, Li Y-X, Ren F-Z, et al. Assessment of the endocrine-disrupting effects of organophosphorus pesticide triazophos and its metabolites on endocrine hormones biosynthesis, transport and receptor binding in Silico. Food and Chemical Toxicology [electronic article]. 2019;133:110759. (https://www.sciencedirect.com/science/article/pii/S0278691519305496). (Accessed February 13, 2024)
54.
Perrier F, Giorgis-Allemand L, Slama R, et al. Within-subject Pooling of Biological Samples to Reduce Exposure Misclassification in Biomarker-based Studies. Epidemiology [electronic article]. 2016;27(3):378. (https://journals.lww.com/epidem/fulltext/2016/05000/within_subject_pooling_of_biological_samples_to.12.aspx). (Accessed March 8, 2024)
55.
Casas M, Basagaña X, Sakhi AK, et al. Variability of urinary concentrations of non-persistent chemicals in pregnant women and school-aged children. Environ. Int. [electronic article]. 2018;121(Pt 1, Pt 1):561–573. (http://dx.doi.org/10.1016/j.envint.2018.09.046)
56.
Rafi Z, Greenland S. Semantic and cognitive tools to aid statistical science: Replace confidence and significance by compatibility and surprise. BMC Medical Research Methodology [electronic article]. 2020;20(1):244. (https://doi.org/10.1186/s12874-020-01105-9). (Accessed March 8, 2024)

Tables for descriptive data

Study populations

Table 1: Participant characteristics (HELIX subcohort; 2013-2016).
Characteristic N = 1,297a
Child age (years) 8.1 (6.5, 8.9)
Child breastfeeding 1,093.0 (84.7%)
    Unknown 6
Child ethnicity
    Caucasian 1,157.0 (90.0%)
    Pakistani 80.0 (6.2%)
    Asian 21.0 (1.6%)
    Other 19.0 (1.5%)
    African 7.0 (0.5%)
    Native American 2.0 (0.2%)
    White non European 0.0 (0.0%)
    Unknown 11
Child head circumference (cm) 51.8 (50.6, 52.9)
    Unknown 3
Child height (m) 1.3 (1.2, 1.4)
Child neuropsychological diagnosis 95.0 (7.3%)
Child rest before assessment
    Yes 1,209.0 (93.3%)
    Not as well as usual 87.0 (6.7%)
    Unknown 1
Child sex
    Male 710.0 (54.7%)
    Female 587.0 (45.3%)
Child weight (kg) 26.9 (22.9, 32.6)
Chiod mood before assessment
    Usual 1,232.0 (95.1%)
    Not usual 64.0 (4.9%)
    Unknown 1
Cohort
    MOBA 272.0 (21.0%)
    INMA 221.0 (17.0%)
    BIB 204.0 (15.7%)
    KANC 203.0 (15.7%)
    RHEA 199.0 (15.3%)
    EDEN 198.0 (15.3%)
Creatinine night sample (g/l) 1.7 (0.9, 3.0)
    Unknown 321
Creatinine pooled sample (g/l) 1.0 (0.8, 1.2)
Date of test (season)
    Spring 358.0 (27.7%)
    Winter 339.0 (26.2%)
    Autumn 300.0 (23.2%)
    Summer 297.0 (23.0%)
    Unknown 3
Family affluence scale
    6 410.0 (31.7%)
    5 325.0 (25.1%)
    7 248.0 (19.2%)
    4 174.0 (13.4%)
    3 92.0 (7.1%)
    2 28.0 (2.2%)
    1 12.0 (0.9%)
    0 6.0 (0.5%)
    Unknown 2
Fast food/take away (times/week) 0.1 (0.1, 0.5)
    Unknown 7
Fasting time before visit (hours) 3.3 (2.8, 4.0)
Financial situation of the parents
    Doing alright 414.0 (32.1%)
    Living comfortably 412.0 (31.9%)
    Getting by 331.0 (25.6%)
    Finding it quite difficult 86.0 (6.7%)
    Finding it very difficult 40.0 (3.1%)
    Does not wish to answer 8.0 (0.6%)
    Unknown 6
Fish and seafood (times/week) 2.0 (1.1, 3.5)
    Unknown 5
Fruits (times/week) 9.0 (5.9, 18.0)
    Unknown 7
Hit reaction time standard error (ms) 299.6 (231.3, 368.2)
    Unknown 18
Marital status
    Living with the father 1,212.0 (94.5%)
    Living alone 39.0 (3.0%)
    Other situation 31.0 (2.4%)
    Unknown 15
Maternal tobacco consumption
    Non-smoker and has never smoked 681.0 (52.6%)
    Daily smoker 200.0 (15.5%)
    Non-smoker but previously smoked daily 186.0 (14.4%)
    Non-smoker but previously smoked although not daily 163.0 (12.6%)
    Smoker but not daily 64.0 (4.9%)
    Unknown 3
Organic food (times/week) 0.5 (0.0, 3.0)
    Unknown 7
Pregnancy maternal active smoking 190.0 (15.1%)
    Unknown 40
Pregnancy maternal passive smoking 514.0 (40.3%)
    Unknown 21
Vegetables (times/week) 6.5 (4.0, 10.0)
    Unknown 6
a n (%); Median (IQR)

Endocrine disruptors

Table 2: Participants endocrine disruptors concentrations expressed in \(\mu\)grams/L (HELIX subcohort; 2013-2016).
Characteristic N = 1,297a N = 1,297b
OP pesticide metabolites
DEP 1.8 (0.4, 4.6) 2.0 (0.2)
DETP 0.1 (0.1, 1.7) 21.0 (1.6)
DMP 0.4 (0.3, 4.6) 6.0 (0.5)
DMTP 2.8 (1.2, 6.3) 1.0 (0.1)
Phenols
BPA 3.8 (2.3, 7.0) 12.0 (0.9)
BUPA 0.1 (0.0, 0.1) 5.0 (0.4)
ETPA 0.7 (0.4, 1.2) 3.0 (0.2)
MEPA 6.3 (3.1, 24.1) 2.0 (0.2)
OXBE 2.0 (0.8, 6.6) 0.0 (0.0)
PRPA 0.2 (0.0, 1.6) 17.0 (1.3)
TRCS 0.6 (0.3, 1.5) 0.0 (0.0)
Phthalate metabolites
MBzP 4.8 (2.7, 8.7) 1.0 (0.1)
MECPP 32.8 (19.9, 57.6) 1.0 (0.1)
MEHHP 19.3 (11.4, 33.1) 3.0 (0.2)
MEHP 2.8 (1.6, 5.1) 41.0 (3.2)
MEOHP 12.2 (7.1, 20.4) 1.0 (0.1)
MEP 32.5 (15.0, 79.2) 0.0 (0.0)
MiBP 40.2 (24.5, 71.1) 0.0 (0.0)
MnBP 22.7 (14.5, 38.8) 0.0 (0.0)
oh-MiNP 5.0 (3.1, 9.3) 0.0 (0.0)
oxo-MiNP 2.7 (1.7, 5.0) 0.0 (0.0)
a Median (IQR)
b N missing (% missing)

Glucocorticosteroids

Table 3: Participants derived glucocorticosteroids concentrations expressed in ng/ml (HELIX subcohort; 2013-2016).
Characteristic N = 1,004a N = 976a,b
cortisol production 4,607.9 (2,860.5, 6,787.6); 18.0 (1.8) 4,559.5 (2,834.5, 6,731.7); 17.0 (1.7)
cortisone production 4,608.1 (2,920.8, 6,843.9); 19.0 (1.9) 4,580.7 (2,899.3, 6,800.5); 18.0 (1.8)
corticosterone production 257.8 (157.9, 410.5); 3.0 (0.3) 256.7 (157.5, 409.7); 3.0 (0.3)
a Median (IQR); N missing (% missing)
b Measurements available for the HELIX subcohort.

Tables for other analyses

Marginal hypotheses for effect modification

Table 4: Pairwise differences between marginal contrasts on the logarithmic scale of males and females, for the effect of a increase from the 10th to the 90th percentile of endocrine disrupting chemicals (EDCs) on hit reaction time standard error (HRT-SE), expressed in ms, and on the glucocorticosteroids, expressed in ng/ml (HELIX subcohort; 2013-2016).
HRT-SEa corticosterone productiona cortisol productiona cortisone productiona
OP pesticide metabolites
DEP 0.019 (-0.022, 0.061) -0.082 (-0.276, 0.113) -0.139 (-0.374, 0.096) -0.104 (-0.312, 0.103)
DETP 0.025 (-0.054, 0.104) -0.16 (-0.332, 0.011) -0.071 (-0.264, 0.123) -0.097 (-0.269, 0.076)
DMP -0.034 (-0.093, 0.025) 0.007 (-0.217, 0.231) -0.031 (-0.119, 0.057) -0.069 (-0.207, 0.07)
DMTP 0.005 (-0.095, 0.106) -0.014 (-0.165, 0.137) -0.21 (-0.326, -0.094) -0.166 (-0.353, 0.022)
Phenols
BPA 0.032 (-0.026, 0.09) -0.153 (-0.291, -0.015) -0.125 (-0.269, 0.018) -0.085 (-0.216, 0.047)
BUPA -0.022 (-0.067, 0.024) -0.117 (-0.247, 0.012) -0.129 (-0.209, -0.048) -0.013 (-0.112, 0.086)
ETPA 0.012 (-0.021, 0.045) -0.254 (-0.416, -0.092) -0.184 (-0.39, 0.022) -0.219 (-0.472, 0.034)
MEPA -0.001 (-0.061, 0.058) -0.129 (-0.271, 0.013) -0.127 (-0.258, 0.004) -0.144 (-0.257, -0.03)
OXBE 0.032 (0.004, 0.061) -0.213 (-0.486, 0.059) -0.077 (-0.306, 0.153) -0.064 (-0.274, 0.146)
PRPA 0.015 (-0.045, 0.074) -0.12 (-0.262, 0.022) -0.043 (-0.238, 0.151) -0.102 (-0.223, 0.019)
TRCS -0.017 (-0.076, 0.042) -0.142 (-0.251, -0.034) -0.13 (-0.248, -0.012) -0.152 (-0.207, -0.096)
Phthalate metabolites
MBzP -0.066 (-0.126, -0.007) -0.026 (-0.098, 0.047) -0.018 (-0.143, 0.108) -0.079 (-0.174, 0.016)
MECPP 0.008 (-0.076, 0.092) -0.014 (-0.165, 0.136) -0.043 (-0.084, -0.002) 0.017 (-0.055, 0.09)
MEHHP 0.028 (-0.075, 0.131) -0.052 (-0.264, 0.161) -0.091 (-0.208, 0.026) -0.006 (-0.087, 0.075)
MEHP 0.017 (-0.082, 0.115) -0.165 (-0.259, -0.071) -0.221 (-0.289, -0.153) -0.177 (-0.298, -0.055)
MEOHP 0.02 (-0.068, 0.107) -0.061 (-0.232, 0.111) -0.075 (-0.157, 0.006) 0.009 (-0.063, 0.08)
MEP -0.053 (-0.138, 0.033) -0.05 (-0.408, 0.308) -0.083 (-0.384, 0.218) -0.119 (-0.338, 0.1)
MiBP -0.02 (-0.138, 0.098) 0.037 (-0.175, 0.249) -0.041 (-0.267, 0.184) -0.021 (-0.162, 0.12)
MnBP -0.035 (-0.11, 0.041) 0.029 (-0.186, 0.243) 0.063 (-0.134, 0.26) 0.017 (-0.076, 0.111)
oh-MiNP 0.046 (-0.009, 0.102) -0.127 (-0.335, 0.08) -0.181 (-0.33, -0.033) -0.164 (-0.304, -0.024)
oxo-MiNP -0.026 (-0.059, 0.008) -0.12 (-0.315, 0.076) -0.146 (-0.303, 0.011) -0.127 (-0.238, -0.016)
a Estimate and 95% CI.
Table 5: Pairwise differences between marginal contrasts on the logarithmic scale of males and females, for the effect of a increase from the 10th to the 90th percentile of the glucocorticosteroids on hit reaction time standard error (HRT-SE) expressed in ms (HELIX subcohort; 2013-2016).
HRT-SEa
Glucocorticosteroids
corticosterone production 0.126 (0.009, 0.243)
cortisol production 0.097 (-0.045, 0.238)
cortisone production 0.14 (0.019, 0.261)
a Estimate and 95% CI.

Figures for main results

Marginal contrasts

Figure 1: Marginal contrasts on the logarithmic scale for the effect of a increase from the 10th to the 90th percentile of the endocrine disrupting chemicals (EDCs) on hit reaction time standard error (HRT-SE) expressed in ms (HELIX subcohort; 2013-2016). Circles indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56).

Figure 2: Marginal contrasts on the logarithmic scale for the effect of a increase from the 10th to the 90th percentile of the endocrine disrupting chemicals (EDCs) on the glucocorticosteroids expressed in ng/ml (HELIX subcohort; 2013-2016). Circles, triangles, and squares indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56).

Figure 3: Marginal contrasts on the logarithmic scale for the effect of a increase from the 10th to the 90th percentile of the glucocorticosteroids on hit reaction time standard error (HRT-SE) expressed in ms (HELIX subcohort; 2013-2016). Circles indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56). Abbreviations: cortisone production (cortisone prod.); cortisol production (cortisol prod.); corticost. prod. (corticosterone production).

Supplementary information

Directed Acyclic Graphs

dag {
age_child
biomarker
breastfeeding
bw
characteristics_child
chemical [exposure]
child_diet
child_smoking
cohort
creatinine
envFactors_visit
ethnicity_child
ethnicity_mother
familySEP
gestational_age
maternalAlcohol_preg
maternalDiet_preg
maternalSEP_preg
maternalSmoking_preg
neuropsychologicalDiagnosis_child
outcome [outcome]
paternalSEP_preg
season_visit
sex_child
time_lastMeal
type_sample
age_child -> biomarker
age_child -> characteristics_child
age_child -> creatinine
age_child -> outcome
age_child -> type_sample
biomarker -> outcome
breastfeeding -> neuropsychologicalDiagnosis_child
breastfeeding -> outcome
bw -> characteristics_child
bw -> neuropsychologicalDiagnosis_child
characteristics_child -> biomarker
characteristics_child -> chemical
characteristics_child -> creatinine
characteristics_child -> outcome
chemical -> biomarker
chemical -> outcome
child_diet -> biomarker
child_diet -> characteristics_child
child_diet -> chemical
child_diet -> outcome
child_smoking -> biomarker
child_smoking -> characteristics_child
child_smoking -> creatinine
child_smoking -> outcome
cohort -> biomarker
cohort -> bw
cohort -> characteristics_child
cohort -> chemical
cohort -> child_diet
cohort -> creatinine
cohort -> outcome
creatinine -> biomarker
creatinine -> chemical
creatinine -> outcome
envFactors_visit -> outcome
ethnicity_child -> biomarker
ethnicity_child -> bw
ethnicity_child -> characteristics_child
ethnicity_child -> chemical
ethnicity_child -> child_diet
ethnicity_child -> child_smoking
ethnicity_child -> creatinine
ethnicity_child -> neuropsychologicalDiagnosis_child
ethnicity_child -> outcome
ethnicity_mother -> biomarker
ethnicity_mother -> breastfeeding
ethnicity_mother -> bw
ethnicity_mother -> characteristics_child
ethnicity_mother -> child_diet
ethnicity_mother -> familySEP
ethnicity_mother -> maternalAlcohol_preg
ethnicity_mother -> maternalDiet_preg
ethnicity_mother -> maternalSEP_preg
ethnicity_mother -> maternalSmoking_preg
ethnicity_mother -> neuropsychologicalDiagnosis_child
ethnicity_mother -> outcome
familySEP -> biomarker
familySEP -> characteristics_child
familySEP -> chemical
familySEP -> child_diet
familySEP -> child_smoking
familySEP -> creatinine
familySEP -> outcome
gestational_age -> bw
gestational_age -> characteristics_child
gestational_age -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> bw
maternalAlcohol_preg -> characteristics_child
maternalAlcohol_preg -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> outcome
maternalDiet_preg -> characteristics_child
maternalDiet_preg -> neuropsychologicalDiagnosis_child
maternalDiet_preg -> outcome
maternalSEP_preg -> breastfeeding
maternalSEP_preg -> bw
maternalSEP_preg -> characteristics_child
maternalSEP_preg -> familySEP
maternalSEP_preg -> maternalAlcohol_preg
maternalSEP_preg -> maternalDiet_preg
maternalSEP_preg -> maternalSmoking_preg
maternalSEP_preg -> neuropsychologicalDiagnosis_child
maternalSEP_preg -> outcome
maternalSmoking_preg -> bw
maternalSmoking_preg -> characteristics_child
maternalSmoking_preg -> neuropsychologicalDiagnosis_child
maternalSmoking_preg -> outcome
neuropsychologicalDiagnosis_child -> outcome
paternalSEP_preg -> breastfeeding
paternalSEP_preg -> bw
paternalSEP_preg -> characteristics_child
paternalSEP_preg -> familySEP
paternalSEP_preg -> maternalAlcohol_preg
paternalSEP_preg -> maternalDiet_preg
paternalSEP_preg -> maternalSmoking_preg
paternalSEP_preg -> neuropsychologicalDiagnosis_child
paternalSEP_preg -> outcome
season_visit -> biomarker
season_visit -> chemical
sex_child -> biomarker
sex_child -> characteristics_child
sex_child -> chemical
sex_child -> child_diet
sex_child -> child_smoking
sex_child -> creatinine
sex_child -> neuropsychologicalDiagnosis_child
sex_child -> outcome
sex_child -> type_sample
time_lastMeal -> biomarker
time_lastMeal -> chemical
type_sample -> chemical
type_sample -> creatinine
}
dag {
age_child
biomarker [outcome]
breastfeeding
bw
characteristics_child
chemical [exposure]
child_diet
child_smoking
cohort
creatinine
envFactors_visit
ethnicity_child
ethnicity_mother
familySEP
gestational_age
maternalAlcohol_preg
maternalDiet_preg
maternalSEP_preg
maternalSmoking_preg
neuropsychologicalDiagnosis_child
outcome
paternalSEP_preg
season_visit
sex_child
time_lastMeal
type_sample
age_child -> biomarker
age_child -> characteristics_child
age_child -> creatinine
age_child -> outcome
age_child -> type_sample
biomarker -> outcome
breastfeeding -> neuropsychologicalDiagnosis_child
breastfeeding -> outcome
bw -> characteristics_child
bw -> neuropsychologicalDiagnosis_child
characteristics_child -> biomarker
characteristics_child -> chemical
characteristics_child -> creatinine
characteristics_child -> outcome
chemical -> biomarker
chemical -> outcome
child_diet -> biomarker
child_diet -> characteristics_child
child_diet -> chemical
child_diet -> outcome
child_smoking -> biomarker
child_smoking -> characteristics_child
child_smoking -> creatinine
child_smoking -> outcome
cohort -> biomarker
cohort -> bw
cohort -> characteristics_child
cohort -> chemical
cohort -> child_diet
cohort -> creatinine
cohort -> outcome
creatinine -> biomarker
creatinine -> chemical
creatinine -> outcome
envFactors_visit -> outcome
ethnicity_child -> biomarker
ethnicity_child -> bw
ethnicity_child -> characteristics_child
ethnicity_child -> chemical
ethnicity_child -> child_diet
ethnicity_child -> child_smoking
ethnicity_child -> creatinine
ethnicity_child -> neuropsychologicalDiagnosis_child
ethnicity_child -> outcome
ethnicity_mother -> biomarker
ethnicity_mother -> breastfeeding
ethnicity_mother -> bw
ethnicity_mother -> characteristics_child
ethnicity_mother -> child_diet
ethnicity_mother -> familySEP
ethnicity_mother -> maternalAlcohol_preg
ethnicity_mother -> maternalDiet_preg
ethnicity_mother -> maternalSEP_preg
ethnicity_mother -> maternalSmoking_preg
ethnicity_mother -> neuropsychologicalDiagnosis_child
ethnicity_mother -> outcome
familySEP -> biomarker
familySEP -> characteristics_child
familySEP -> chemical
familySEP -> child_diet
familySEP -> child_smoking
familySEP -> creatinine
familySEP -> outcome
gestational_age -> bw
gestational_age -> characteristics_child
gestational_age -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> bw
maternalAlcohol_preg -> characteristics_child
maternalAlcohol_preg -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> outcome
maternalDiet_preg -> characteristics_child
maternalDiet_preg -> neuropsychologicalDiagnosis_child
maternalDiet_preg -> outcome
maternalSEP_preg -> breastfeeding
maternalSEP_preg -> bw
maternalSEP_preg -> characteristics_child
maternalSEP_preg -> familySEP
maternalSEP_preg -> maternalAlcohol_preg
maternalSEP_preg -> maternalDiet_preg
maternalSEP_preg -> maternalSmoking_preg
maternalSEP_preg -> neuropsychologicalDiagnosis_child
maternalSEP_preg -> outcome
maternalSmoking_preg -> bw
maternalSmoking_preg -> characteristics_child
maternalSmoking_preg -> neuropsychologicalDiagnosis_child
maternalSmoking_preg -> outcome
neuropsychologicalDiagnosis_child -> outcome
paternalSEP_preg -> breastfeeding
paternalSEP_preg -> bw
paternalSEP_preg -> characteristics_child
paternalSEP_preg -> familySEP
paternalSEP_preg -> maternalAlcohol_preg
paternalSEP_preg -> maternalDiet_preg
paternalSEP_preg -> maternalSmoking_preg
paternalSEP_preg -> neuropsychologicalDiagnosis_child
paternalSEP_preg -> outcome
season_visit -> biomarker
season_visit -> chemical
sex_child -> biomarker
sex_child -> characteristics_child
sex_child -> chemical
sex_child -> child_diet
sex_child -> child_smoking
sex_child -> creatinine
sex_child -> neuropsychologicalDiagnosis_child
sex_child -> outcome
sex_child -> type_sample
time_lastMeal -> biomarker
time_lastMeal -> chemical
type_sample -> chemical
type_sample -> creatinine
}
dag {
age_child
biomarker [exposure]
breastfeeding
bw
characteristics_child
chemical
child_diet
child_smoking
cohort
creatinine
envFactors_visit
ethnicity_child
ethnicity_mother
familySEP
gestational_age
maternalAlcohol_preg
maternalDiet_preg
maternalSEP_preg
maternalSmoking_preg
neuropsychologicalDiagnosis_child
outcome [outcome]
paternalSEP_preg
season_visit
sex_child
time_lastMeal
type_sample
age_child -> biomarker
age_child -> characteristics_child
age_child -> creatinine
age_child -> outcome
age_child -> type_sample
biomarker -> outcome
breastfeeding -> neuropsychologicalDiagnosis_child
breastfeeding -> outcome
bw -> characteristics_child
bw -> neuropsychologicalDiagnosis_child
characteristics_child -> biomarker
characteristics_child -> chemical
characteristics_child -> creatinine
characteristics_child -> outcome
chemical -> biomarker
chemical -> outcome
child_diet -> biomarker
child_diet -> characteristics_child
child_diet -> chemical
child_diet -> outcome
child_smoking -> biomarker
child_smoking -> characteristics_child
child_smoking -> creatinine
child_smoking -> outcome
cohort -> biomarker
cohort -> bw
cohort -> characteristics_child
cohort -> chemical
cohort -> child_diet
cohort -> creatinine
cohort -> outcome
creatinine -> biomarker
creatinine -> chemical
creatinine -> outcome
envFactors_visit -> outcome
ethnicity_child -> biomarker
ethnicity_child -> bw
ethnicity_child -> characteristics_child
ethnicity_child -> chemical
ethnicity_child -> child_diet
ethnicity_child -> child_smoking
ethnicity_child -> creatinine
ethnicity_child -> neuropsychologicalDiagnosis_child
ethnicity_child -> outcome
ethnicity_mother -> biomarker
ethnicity_mother -> breastfeeding
ethnicity_mother -> bw
ethnicity_mother -> characteristics_child
ethnicity_mother -> child_diet
ethnicity_mother -> familySEP
ethnicity_mother -> maternalAlcohol_preg
ethnicity_mother -> maternalDiet_preg
ethnicity_mother -> maternalSEP_preg
ethnicity_mother -> maternalSmoking_preg
ethnicity_mother -> neuropsychologicalDiagnosis_child
ethnicity_mother -> outcome
familySEP -> biomarker
familySEP -> characteristics_child
familySEP -> chemical
familySEP -> child_diet
familySEP -> child_smoking
familySEP -> creatinine
familySEP -> outcome
gestational_age -> bw
gestational_age -> characteristics_child
gestational_age -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> bw
maternalAlcohol_preg -> characteristics_child
maternalAlcohol_preg -> neuropsychologicalDiagnosis_child
maternalAlcohol_preg -> outcome
maternalDiet_preg -> characteristics_child
maternalDiet_preg -> neuropsychologicalDiagnosis_child
maternalDiet_preg -> outcome
maternalSEP_preg -> breastfeeding
maternalSEP_preg -> bw
maternalSEP_preg -> characteristics_child
maternalSEP_preg -> familySEP
maternalSEP_preg -> maternalAlcohol_preg
maternalSEP_preg -> maternalDiet_preg
maternalSEP_preg -> maternalSmoking_preg
maternalSEP_preg -> neuropsychologicalDiagnosis_child
maternalSEP_preg -> outcome
maternalSmoking_preg -> bw
maternalSmoking_preg -> characteristics_child
maternalSmoking_preg -> neuropsychologicalDiagnosis_child
maternalSmoking_preg -> outcome
neuropsychologicalDiagnosis_child -> outcome
paternalSEP_preg -> breastfeeding
paternalSEP_preg -> bw
paternalSEP_preg -> characteristics_child
paternalSEP_preg -> familySEP
paternalSEP_preg -> maternalAlcohol_preg
paternalSEP_preg -> maternalDiet_preg
paternalSEP_preg -> maternalSmoking_preg
paternalSEP_preg -> neuropsychologicalDiagnosis_child
paternalSEP_preg -> outcome
season_visit -> biomarker
season_visit -> chemical
sex_child -> biomarker
sex_child -> characteristics_child
sex_child -> chemical
sex_child -> child_diet
sex_child -> child_smoking
sex_child -> creatinine
sex_child -> neuropsychologicalDiagnosis_child
sex_child -> outcome
sex_child -> type_sample
time_lastMeal -> biomarker
time_lastMeal -> chemical
type_sample -> chemical
type_sample -> creatinine
}

Supplementary tables

Tables for descriptive data

Information about the endocrine disruptors

Compound Symbol Variable name PubChem CID Parental compound
OP pesticide metabolites
diethyl dithiophosphate DEDTP dedtp 9274
diethyl phosphate DEP dep 654
diethyl thiophosphate DETP detp 3683036
dimethyl dithiophosphate DMDTP dmdtp
dimethyl phosphate DMP dmp 13134
dimethyl thiophosphate DMTP dmtp 168140
Phenols
bisphenol A BPA bpa 6623
n‑butyl‑paraben BUPA bupa 7184
ethyl-paraben ETPA etpa 8434
methyl-paraben MEPA mepa 7456
oxybenzone OXBE oxbe 4632
propyl-paraben PRPA prpa 7175
triclosan TRCS trcs 5564
Phthalate metabolites
mono benzyl phthalate MBzP mbzp 31736 BBzP
mono‑2‑ethyl 5‑carboxypentyl phthalate MECPP mecpp 148386 DEHP
mono‑2‑ethyl‑5‑hydroxyhexyl phthalate MEHHP mehhp 170295 DEHP
mono‑2‑ethylhexyl phthalate MEHP mehp 21924291 DEHP
mono‑2‑ethyl‑5‑oxohexyl phthalate MEOHP meohp 119096 DEHP
monoethyl phthalate MEP mep 75318 DEP
mono‑iso‑butyl phthalate MiBP mibp 92272 DiBP
mono‑n‑butyl phthalate MnBP mnbp 8575 DnBP
mono‑4‑methyl‑7‑hydroxyoctyl phthalate oh-MiNP ohminp 102401880 MiNP
mono‑4‑methyl‑7‑oxooctyl phthalate oxo-MiNP oxominp 102401881 MiNP
Table S1: Information about non-persistent endocrine disrupting chemicals (EDCs), including the full compound name, the standard symbol, the used variable name, the identifier from PubChem, and the parental compound.

Information about the glucocorticosteroids

Metabolite Symbol HMDB ID CAS number
Androgen
Androsternedione AED HMDB0000053 63-05-8
Testosterone T HMDB0000234 58-22-0
Androgen metabolite
Androsterone Andros HMDB0000031 53-41-8
Etiocholanolone Etio HMDB0000490 53-42-9
Glucocorticosteroid
11-dehydrocorticosterone A HMDB0004029 72-23-1
Corticosterone B HMDB0001547 50-22-6
Cortisol F HMDB0000063 50-23-7
Cortisone E HMDB0002802 53-06-5
Glucocorticosteroid metabolite
11β-hydroxyandrosterone 11OHAndros HMDB0002984 57-61-4
17-deoxycortolone 17-DO-cortolone NA NA
20α-dihydrocortisol 20aDHF NA NA
20α-dihydrocortisone 20aDHE NA NA
20β-dihydrocortisol 20bDHF NA NA
20β-dihydrocortisone 20bDHE NA NA
5α,20α-cortol 5a20acortol HMDB0003180 516-38-1
5α,20β-cortol 5a20bcortol HMDB0005821 667-65-2
5α-tetrahydrocorticosterone 5aTHB HMDB0000449 600-63-5
5α-tetrahydrocortisol 5aTHF HMDB0000526 302-91-0
5α-tetrahydrocortisone 5aTHE NA NA
5β,20α-cortol 5b20acortol HMDB0003180 516-38-1
5β,20α-cortolone 5b20acortolone HMDB0003128 516-42-7
5β,20β-cortol 5b20bcortol HMDB0005821 667-65-2
5β,20β-cortolone 5b20bcortolone NA NA
5β-dihydrocortisol 5bDHF HMDB0003259 1482-50-4
5β-tetrahydrocorticosterone 5bTHB HMDB0000268 68-42-8
5β-tetrahydrocortisol 5bTHF HMDB0000949 1953-02-01
5β-tetrahydrocortisone 5bTHE NA NA
6β-hydroxycortisol 6OHF HMDB0247074
6β-hydroxycortisone 6OHE NA NA
Glucocorticosteroid precursor
17-hydroxyprogesterone 17OHP HMDB0000374 68-96-2
Cortexolone S HMDB0000015 152-58-9
Deoxycorticosterone DOC HMDB0000016 64-85-7
Glucocorticosteroid precursor metabolite
17-hydroxypregnanolone 17HP HMDB0000363 387-79-1
5β-dihydrocortexolone 5bDHS NA NA
5β-tetrahydrocortexolone 5bTHS NA NA
Pregnantriol PT NA 1098-45-9
Tetrahydrocortexolone THS HMDB0005972 68-60-0
Abbreviations: Human Metabolome Database (HMDB); Chemical Abstracts Service (CAS).
Table S2: Information about the glucocorticosteroids, including the full metabolite name, the standard symbol, the identifier from the HMDB, and the CAS number.

Codebooks

type description coding labels remarks comments includeda
age_child
hs_age_years numerical Age years TRUE
breastfeeding
hs_bf categorical Child breastfeeding 0,1 No, Yes TRUE
characteristics_child
hs_c_height numerical Height m TRUE
hs_c_weight numerical Weight kg TRUE
hs_head_circ numerical Head circumference cm TRUE
child_diet
hs_fastfood numerical Fast food/take away Times / week TRUE
hs_org_food numerical Organic food Times / week TRUE
hs_total_fish numerical Fish and seafood Times / week TRUE
hs_total_fruits numerical Fruits Times / week TRUE
hs_total_veg numerical Vegetables Times / week TRUE
child_smoking
hs_tob categorical Tobacco consumption 1,2,3,4,5 Non-smoker and has never smoked, Non-smoker but previously smoked although not daily, Non-smoker but previously smoked daily, Smoker but not daily, Daily smoker TRUE
cohort
cohort character Cohort SAB,EDEN,BIB,RHEA,KANC,MOBA SAB, EDEN, BIB, RHEA, KANC, MOBA TRUE
creatinine
hs_creatinine_cg numerical Creatinine pooled sample Values below the limit of detection imputed G / L TRUE
envFactors_visit
hs_mood categorical Mood before assessment 1,2 Usual, Not usual TRUE
hs_rest_nth categorical Rest before assessment 1,2 Yes, Not as well as usual TRUE
ethnicity_child
h_ethnicity_c character Child ethnicity 1,2,3,4,5,6,7 African, Asian, Caucasian, Native American, Other, Pakistani, White non European TRUE
ethnicity_mother
h_ethnicity_m integer Mother ethnicity 1,2,3,4,5,6,7 White European, Pakistani, Asian, African, Other, Native American, White non European FALSE
familySEP
FAS_score numerical Family Affluence Scale TRUE
hs_finance categorical Financial situation 1,2,3,4,5,6 Living comfortably, Doing alright, Getting by, Finding it quite difficult, Finding it very difficult, Does not wish to answer TRUE
maternalAlcohol_preg
e3_alcpreg_g numerical Alcool during pregnancy Glasses / week FALSE
maternalDiet_preg
h_cereal_preg numerical Cereal consumption during pregnancy Times / week FALSE
h_dairy_preg numerical Dairy consumption during pregnancy Times / week FALSE
h_fastfood_preg numerical Fast food consumption during pregnancy Times / week FALSE
h_fish_preg numerical Fish consumption during pregnancy Times / week FALSE
h_fruit_preg numerical Fruit consumption during pregnancy Times / week FALSE
h_legume_preg numerical Legume consumption during pregnancy Times / week FALSE
h_meat_preg numerical Meat consumption during pregnancy Times / week FALSE
h_veg_preg numerical Vegetables consumption during pregnancy Times / week FALSE
maternalSEP_preg
e3_edum categorical Maternal education 0,1,2 Primary school, Secondary school, University degree or higher FALSE
e3_marital categorical Marital status 0,1,2 Living with the father, Living alone, Other situation TRUE
e3_ses categorical Socioeconomic status of the parents 1,2,3 Low income, Medium income, High income FALSE
maternalSmoking_preg
e3_asmokyn_p categorical Pregnancy maternal active smoking 0,1 No, Yes TRUE
e3_psmokanyt categorical Pregnancy maternal passive smoking 0,1 No, Yes TRUE
neuropsychologicalDiagnosis_child
hs_neuro_diag categorical Child neuropsychological diagnosis 1,2 No, Yes TRUE
paternalSEP_preg
e3_eduf categorical Paternal education 0,1,2 Primary school, Secondary school, University degree or higher FALSE
season_visit
hs_date_neu date Date of test season TRUE
sex_child
e3_sex categorical Sex 0,1 Male, Female TRUE
time_lastMeal
hs_dift_mealblood_imp numerical Fasting time hours TRUE
a Percentage of confounders included in the models: 65.79%.
Table S3: Codebook for the covariates used in the estimation of the marginal comparisons of endocrine disrupting chemicals (EDCs) on hit reaction time standard error (HRT-SE).
type description coding labels remarks comments includeda
age_child
hs_age_years numerical Age years TRUE
characteristics_child
hs_c_height numerical Height m TRUE
hs_c_weight numerical Weight kg TRUE
hs_head_circ numerical Head circumference cm TRUE
child_diet
hs_fastfood numerical Fast food/take away Times / week TRUE
hs_org_food numerical Organic food Times / week TRUE
hs_total_fish numerical Fish and seafood Times / week TRUE
hs_total_fruits numerical Fruits Times / week TRUE
hs_total_veg numerical Vegetables Times / week TRUE
child_smoking
hs_tob categorical Tobacco consumption 1,2,3,4,5 Non-smoker and has never smoked, Non-smoker but previously smoked although not daily, Non-smoker but previously smoked daily, Smoker but not daily, Daily smoker TRUE
cohort
cohort character Cohort SAB,EDEN,BIB,RHEA,KANC,MOBA SAB, EDEN, BIB, RHEA, KANC, MOBA TRUE
creatinine
creatinine_to_helix numerical Creatinine night sample G / L TRUE
hs_creatinine_cg numerical Creatinine pooled sample Values below the limit of detection imputed G / L TRUE
ethnicity_child
h_ethnicity_c character Child ethnicity 1,2,3,4,5,6,7 African, Asian, Caucasian, Native American, Other, Pakistani, White non European TRUE
ethnicity_mother
h_ethnicity_m integer Mother ethnicity 1,2,3,4,5,6,7 White European, Pakistani, Asian, African, Other, Native American, White non European FALSE
familySEP
FAS_score numerical Family Affluence Scale TRUE
hs_finance categorical Financial situation 1,2,3,4,5,6 Living comfortably, Doing alright, Getting by, Finding it quite difficult, Finding it very difficult, Does not wish to answer TRUE
season_visit
hs_date_neu date Date of test season TRUE
sex_child
e3_sex categorical Sex 0,1 Male, Female TRUE
time_lastMeal
hs_dift_mealblood_imp numerical Fasting time hours TRUE
a Percentage of confounders included in the models: 95%.
Table S4: Codebook for the covariates used in the estimation of the marginal comparisons of endocrine disrupting chemicals (EDCs) on the glucocorticosteroids.
type description coding labels remarks comments includeda
age_child
hs_age_years numerical Age years TRUE
breastfeeding
hs_bf categorical Child breastfeeding 0,1 No, Yes TRUE
characteristics_child
hs_c_height numerical Height m TRUE
hs_c_weight numerical Weight kg TRUE
hs_head_circ numerical Head circumference cm TRUE
chemical
hs_bpa_c numerical Bisphenol A (BPA) Values below the limit of detection imputed microg / L TRUE
hs_bupa_c numerical N-Butyl paraben (BUPA) Values below the limit of detection imputed microg / L TRUE
hs_dedtp_cadj numerical Diethyl dithiophosphate (DEDTP) adjusted for creatinine Values below the limit of detection imputed microg / g FALSE
hs_dep_c numerical Diethyl phosphate (DEP) Values below the limit of detection imputed microg / L TRUE
hs_detp_c numerical Diethyl thiophosphate (DETP) Values below the limit of detection imputed microg / L TRUE
hs_dmdtp_craw numerical Dimethyl dithiophosphate (DMDTP) Values below the limit of detection imputed microg / L FALSE
hs_dmp_c numerical Dimethyl phosphate (DMP) Values below the limit of detection imputed microg / L TRUE
hs_dmtp_c numerical Dimethyl thiophosphate (DMTP) Values below the limit of detection imputed microg / L TRUE
hs_etpa_c numerical Ethyl paraben (ETPA) Values below the limit of detection imputed microg / L TRUE
hs_mbzp_c numerical Mono benzyl phthalate (MbzP) Values below the limit of detection imputed microg / L TRUE
hs_mecpp_c numerical Mono-2-ethyl 5-carboxypentyl phthalate (MECPP) Values below the limit of detection imputed microg / L TRUE
hs_mehhp_c numerical Mono-2-ethyl-5-hydroxyhexyl phthalate (MEHHP) Values below the limit of detection imputed microg / L TRUE
hs_mehp_c numerical Mono-2-ethylhexyl phthalate (MEHP) Values below the limit of detection imputed microg / L TRUE
hs_meohp_c numerical Mono-2-ethyl-5-oxohexyl phthalate (MEOHP) Values below the limit of detection imputed microg / L TRUE
hs_mep_c numerical Monoethyl phthalate (MEP) Values below the limit of detection imputed microg / L TRUE
hs_mepa_c numerical Methyl paraben (MEPA) Values below the limit of detection imputed microg / L TRUE
hs_mibp_c numerical Mono-iso-butyl phthalate (MiBP) Values below the limit of detection imputed microg / L TRUE
hs_mnbp_c numerical Mono-n-butyl phthalate (MnBP) Values below the limit of detection imputed microg / L TRUE
hs_ohminp_c numerical Mono-4-methyl-7-hydroxyoctyl phthalate (OHMiNP) Values below the limit of detection imputed microg / L TRUE
hs_oxbe_c numerical Oxybenzone (OXBE) Values below the limit of detection imputed microg / L TRUE
hs_oxominp_c numerical Mono-4-methyl-7-oxooctyl phthalate (OXOMiNP) Values below the limit of detection imputed microg / L TRUE
hs_prpa_c numerical Propyl paraben (PRPA) Values below the limit of detection imputed microg / L TRUE
hs_trcs_c numerical Triclosan (TRCS) Values below the limit of detection imputed microg / L TRUE
child_diet
hs_fastfood numerical Fast food/take away Times / week TRUE
hs_org_food numerical Organic food Times / week TRUE
hs_total_fish numerical Fish and seafood Times / week TRUE
hs_total_fruits numerical Fruits Times / week TRUE
hs_total_veg numerical Vegetables Times / week TRUE
child_smoking
hs_tob categorical Tobacco consumption 1,2,3,4,5 Non-smoker and has never smoked, Non-smoker but previously smoked although not daily, Non-smoker but previously smoked daily, Smoker but not daily, Daily smoker TRUE
cohort
cohort character Cohort SAB,EDEN,BIB,RHEA,KANC,MOBA SAB, EDEN, BIB, RHEA, KANC, MOBA TRUE
creatinine
creatinine_to_helix numerical Creatinine night sample G / L TRUE
envFactors_visit
hs_mood categorical Mood before assessment 1,2 Usual, Not usual TRUE
hs_rest_nth categorical Rest before assessment 1,2 Yes, Not as well as usual TRUE
ethnicity_child
h_ethnicity_c character Child ethnicity 1,2,3,4,5,6,7 African, Asian, Caucasian, Native American, Other, Pakistani, White non European TRUE
ethnicity_mother
h_ethnicity_m integer Mother ethnicity 1,2,3,4,5,6,7 White European, Pakistani, Asian, African, Other, Native American, White non European FALSE
familySEP
FAS_score numerical Family Affluence Scale TRUE
hs_finance categorical Financial situation 1,2,3,4,5,6 Living comfortably, Doing alright, Getting by, Finding it quite difficult, Finding it very difficult, Does not wish to answer TRUE
maternalAlcohol_preg
e3_alcpreg_g numerical Alcool during pregnancy Glasses / week FALSE
maternalDiet_preg
h_cereal_preg numerical Cereal consumption during pregnancy Times / week FALSE
h_dairy_preg numerical Dairy consumption during pregnancy Times / week FALSE
h_fastfood_preg numerical Fast food consumption during pregnancy Times / week FALSE
h_fish_preg numerical Fish consumption during pregnancy Times / week FALSE
h_fruit_preg numerical Fruit consumption during pregnancy Times / week FALSE
h_legume_preg numerical Legume consumption during pregnancy Times / week FALSE
h_meat_preg numerical Meat consumption during pregnancy Times / week FALSE
h_veg_preg numerical Vegetables consumption during pregnancy Times / week FALSE
maternalSEP_preg
e3_edum categorical Maternal education 0,1,2 Primary school, Secondary school, University degree or higher FALSE
e3_marital categorical Marital status 0,1,2 Living with the father, Living alone, Other situation TRUE
e3_ses categorical Socioeconomic status of the parents 1,2,3 Low income, Medium income, High income FALSE
maternalSmoking_preg
e3_asmokyn_p categorical Pregnancy maternal active smoking 0,1 No, Yes TRUE
e3_psmokanyt categorical Pregnancy maternal passive smoking 0,1 No, Yes TRUE
neuropsychologicalDiagnosis_child
hs_neuro_diag categorical Child neuropsychological diagnosis 1,2 No, Yes TRUE
paternalSEP_preg
e3_eduf categorical Paternal education 0,1,2 Primary school, Secondary school, University degree or higher FALSE
sex_child
e3_sex categorical Sex 0,1 Male, Female TRUE
a Percentage of confounders included in the models: 74.58%.
Table S5: Codebook for the covariates used in the estimation of the marginal comparisons of the glucocorticosteroids on hit reaction time standard error (HRT-SE).

Lower limits of quantification of the glucocorticosteroids

Metabolite LLOQ
5aTHF 5.00
5bTHE 5.00
5b20acortolone 5.00
5b20bcortolone 5.00
5a20acortol 2.50
5a20bcortol 2.50
5b20acortol 2.50
5b20bcortol 2.50
11OHAndros 2.00
17HP 2.00
PT 2.00
20bDHF 0.50
5bTHF 0.50
6OHF 0.50
E 0.50
20aDHE 0.50
20bDHE 0.50
5aTHE 0.50
6OHE 0.50
5aTHB 0.50
5bTHB 0.50
17DOcortolone 0.50
5bTHS 0.50
Andros 0.50
Etio 0.50
F 0.25
20aDHF 0.25
5bDHF 0.10
A 0.10
S 0.10
5bDHS 0.10
T 0.10
AED 0.10
Abbreviations: lower limit of quantification (LLOQ).
Table S6: Lower limits of quantification expressed in ng/ml for the glucocorticosteroids (HELIX subcohort; 2013-2016).

Study populations

Characteristic Overall, N = 1,297a BIB, N = 204a EDEN, N = 198a INMA, N = 221a KANC, N = 203a MOBA, N = 272a RHEA, N = 199a
Child age (years) 8.1 (6.5, 8.9) 6.6 (6.5, 6.8) 10.9 (10.4, 11.2) 8.8 (8.4, 9.3) 6.4 (6.1, 6.9) 8.5 (8.2, 8.8) 6.5 (6.4, 6.6)
Child breastfeeding 1,093.0 (84.7%) 147.0 (72.4%) 128.0 (65.0%) 195.0 (88.6%) 187.0 (92.6%) 260.0 (96.3%) 176.0 (88.4%)
    Unknown 6 1 1 1 1 2 0
Child ethnicity






    Caucasian 1,157.0 (90.0%) 87.0 (42.6%) 196.0 (99.5%) 221.0 (100.0%) 200.0 (100.0%) 254.0 (95.8%) 199.0 (100.0%)
    Pakistani 80.0 (6.2%) 80.0 (39.2%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%)
    Asian 21.0 (1.6%) 13.0 (6.4%) 1.0 (0.5%) 0.0 (0.0%) 0.0 (0.0%) 7.0 (2.6%) 0.0 (0.0%)
    Other 19.0 (1.5%) 17.0 (8.3%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 2.0 (0.8%) 0.0 (0.0%)
    African 7.0 (0.5%) 7.0 (3.4%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%)
    Native American 2.0 (0.2%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 2.0 (0.8%) 0.0 (0.0%)
    White non European 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%)
    Unknown 11 0 1 0 3 7 0
Child head circumference (cm) 51.8 (50.6, 52.9) 51.4 (50.3, 52.3) 50.5 (49.5, 52.0) 52.3 (51.3, 53.3) 52.0 (51.0, 53.0) 52.5 (51.5, 53.6) 51.2 (50.2, 52.0)
    Unknown 3 0 0 0 0 0 3
Child height (m) 1.3 (1.2, 1.4) 1.2 (1.2, 1.2) 1.4 (1.4, 1.5) 1.3 (1.3, 1.4) 1.2 (1.2, 1.3) 1.3 (1.3, 1.4) 1.2 (1.2, 1.2)
Child neuropsychological diagnosis 95.0 (7.3%) 3.0 (1.5%) 58.0 (29.3%) 24.0 (10.9%) 1.0 (0.5%) 1.0 (0.4%) 8.0 (4.0%)
Child rest before assessment






    Yes 1,209.0 (93.3%) 192.0 (94.1%) 170.0 (86.3%) 206.0 (93.2%) 200.0 (98.5%) 259.0 (95.2%) 182.0 (91.5%)
    Not as well as usual 87.0 (6.7%) 12.0 (5.9%) 27.0 (13.7%) 15.0 (6.8%) 3.0 (1.5%) 13.0 (4.8%) 17.0 (8.5%)
    Unknown 1 0 1 0 0 0 0
Child sex






    Male 710.0 (54.7%) 112.0 (54.9%) 113.0 (57.1%) 120.0 (54.3%) 111.0 (54.7%) 143.0 (52.6%) 111.0 (55.8%)
    Female 587.0 (45.3%) 92.0 (45.1%) 85.0 (42.9%) 101.0 (45.7%) 92.0 (45.3%) 129.0 (47.4%) 88.0 (44.2%)
Child weight (kg) 26.9 (22.9, 32.6) 22.3 (20.3, 25.0) 35.7 (32.4, 41.2) 30.7 (26.8, 36.5) 23.6 (21.4, 27.1) 28.5 (25.7, 31.6) 23.3 (21.2, 27.2)
Chiod mood before assessment






    Usual 1,232.0 (95.1%) 198.0 (97.1%) 176.0 (89.3%) 214.0 (96.8%) 187.0 (92.1%) 262.0 (96.3%) 195.0 (98.0%)
    Not usual 64.0 (4.9%) 6.0 (2.9%) 21.0 (10.7%) 7.0 (3.2%) 16.0 (7.9%) 10.0 (3.7%) 4.0 (2.0%)
    Unknown 1 0 1 0 0 0 0
Creatinine night sample (g/l) 1.7 (0.9, 3.0) 0.8 (0.6, 1.1) 3.3 (2.0, 4.3) 2.5 (1.5, 3.8) 1.7 (0.9, 2.7) 2.0 (1.2, 3.0) 0.8 (0.4, 1.3)
    Unknown 321 72 64 19 23 72 71
Creatinine pooled sample (g/l) 1.0 (0.8, 1.2) 1.0 (0.8, 1.2) 1.2 (1.0, 1.5) 1.0 (0.8, 1.3) 0.9 (0.7, 1.1) 0.9 (0.7, 1.1) 0.9 (0.7, 1.1)
Date of test (season)






    Spring 358.0 (27.7%) 48.0 (23.5%) 64.0 (32.3%) 71.0 (32.4%) 61.0 (30.0%) 37.0 (13.6%) 77.0 (38.9%)
    Winter 339.0 (26.2%) 40.0 (19.6%) 61.0 (30.8%) 97.0 (44.3%) 38.0 (18.7%) 73.0 (26.8%) 30.0 (15.2%)
    Autumn 300.0 (23.2%) 49.0 (24.0%) 1.0 (0.5%) 30.0 (13.7%) 77.0 (37.9%) 105.0 (38.6%) 38.0 (19.2%)
    Summer 297.0 (23.0%) 67.0 (32.8%) 72.0 (36.4%) 21.0 (9.6%) 27.0 (13.3%) 57.0 (21.0%) 53.0 (26.8%)
    Unknown 3 0 0 2 0 0 1
Family affluence scale






    6 410.0 (31.7%) 34.0 (16.7%) 64.0 (32.3%) 75.0 (34.1%) 50.0 (24.8%) 142.0 (52.2%) 45.0 (22.6%)
    5 325.0 (25.1%) 48.0 (23.5%) 29.0 (14.6%) 65.0 (29.5%) 69.0 (34.2%) 57.0 (21.0%) 57.0 (28.6%)
    7 248.0 (19.2%) 26.0 (12.7%) 90.0 (45.5%) 43.0 (19.5%) 14.0 (6.9%) 53.0 (19.5%) 22.0 (11.1%)
    4 174.0 (13.4%) 40.0 (19.6%) 13.0 (6.6%) 22.0 (10.0%) 38.0 (18.8%) 16.0 (5.9%) 45.0 (22.6%)
    3 92.0 (7.1%) 34.0 (16.7%) 2.0 (1.0%) 11.0 (5.0%) 22.0 (10.9%) 3.0 (1.1%) 20.0 (10.1%)
    2 28.0 (2.2%) 16.0 (7.8%) 0.0 (0.0%) 1.0 (0.5%) 4.0 (2.0%) 0.0 (0.0%) 7.0 (3.5%)
    1 12.0 (0.9%) 4.0 (2.0%) 0.0 (0.0%) 2.0 (0.9%) 4.0 (2.0%) 1.0 (0.4%) 1.0 (0.5%)
    0 6.0 (0.5%) 2.0 (1.0%) 0.0 (0.0%) 1.0 (0.5%) 1.0 (0.5%) 0.0 (0.0%) 2.0 (1.0%)
    Unknown 2 0 0 1 1 0 0
Fast food/take away (times/week) 0.1 (0.1, 0.5) 0.5 (0.1, 1.0) 0.1 (0.1, 0.5) 0.1 (0.1, 0.5) 0.1 (0.0, 0.1) 0.1 (0.1, 0.5) 0.5 (0.1, 0.5)
    Unknown 7 0 0 5 2 0 0
Fasting time before visit (hours) 3.3 (2.8, 4.0) 3.3 (2.8, 4.1) 3.2 (2.8, 3.7) 3.0 (2.6, 3.8) 3.3 (2.8, 3.8) 3.4 (2.8, 3.8) 4.0 (3.3, 4.8)
Financial situation of the parents






    Doing alright 414.0 (32.1%) 73.0 (35.8%) 94.0 (47.5%) 64.0 (29.2%) 61.0 (30.5%) 64.0 (23.5%) 58.0 (29.3%)
    Living comfortably 412.0 (31.9%) 59.0 (28.9%) 49.0 (24.7%) 29.0 (13.2%) 48.0 (24.0%) 202.0 (74.3%) 25.0 (12.6%)
    Getting by 331.0 (25.6%) 59.0 (28.9%) 36.0 (18.2%) 82.0 (37.4%) 70.0 (35.0%) 4.0 (1.5%) 80.0 (40.4%)
    Finding it quite difficult 86.0 (6.7%) 8.0 (3.9%) 9.0 (4.5%) 29.0 (13.2%) 12.0 (6.0%) 1.0 (0.4%) 27.0 (13.6%)
    Finding it very difficult 40.0 (3.1%) 5.0 (2.5%) 10.0 (5.1%) 15.0 (6.8%) 2.0 (1.0%) 0.0 (0.0%) 8.0 (4.0%)
    Does not wish to answer 8.0 (0.6%) 0.0 (0.0%) 0.0 (0.0%) 0.0 (0.0%) 7.0 (3.5%) 1.0 (0.4%) 0.0 (0.0%)
    Unknown 6 0 0 2 3 0 1
Fish and seafood (times/week) 2.0 (1.1, 3.5) 2.0 (1.0, 3.1) 2.1 (1.4, 3.0) 3.5 (2.1, 5.0) 1.0 (0.4, 1.6) 2.6 (1.6, 5.0) 1.5 (1.0, 2.0)
    Unknown 5 1 0 2 2 0 0
Fruits (times/week) 9.0 (5.9, 18.0) 15.5 (10.0, 21.0) 6.6 (3.3, 13.5) 7.5 (3.6, 12.6) 7.3 (3.8, 9.6) 14.1 (8.6, 21.0) 8.5 (6.2, 13.5)
    Unknown 7 2 0 2 2 1 0
Hit reaction time standard error (ms) 299.6 (231.3, 368.2) 355.1 (292.1, 397.5) 237.7 (184.7, 307.0) 256.0 (197.4, 313.8) 368.4 (324.2, 406.6) 248.7 (193.0, 300.9) 340.9 (281.1, 399.2)
    Unknown 18 3 11 3 0 0 1
Marital status






    Living with the father 1,212.0 (94.5%) 178.0 (87.3%) 193.0 (98.0%) 219.0 (99.1%) 168.0 (84.4%) 260.0 (98.5%) 194.0 (98.5%)
    Living alone 39.0 (3.0%) 0.0 (0.0%) 2.0 (1.0%) 0.0 (0.0%) 31.0 (15.6%) 3.0 (1.1%) 3.0 (1.5%)
    Other situation 31.0 (2.4%) 26.0 (12.7%) 2.0 (1.0%) 2.0 (0.9%) 0.0 (0.0%) 1.0 (0.4%) 0.0 (0.0%)
    Unknown 15 0 1 0 4 8 2
Maternal tobacco consumption






    Non-smoker and has never smoked 681.0 (52.6%) 148.0 (72.5%) 87.0 (43.9%) 103.0 (46.8%) 104.0 (51.7%) 138.0 (50.7%) 101.0 (50.8%)
    Daily smoker 200.0 (15.5%) 27.0 (13.2%) 45.0 (22.7%) 45.0 (20.5%) 24.0 (11.9%) 6.0 (2.2%) 53.0 (26.6%)
    Non-smoker but previously smoked daily 186.0 (14.4%) 11.0 (5.4%) 37.0 (18.7%) 42.0 (19.1%) 21.0 (10.4%) 53.0 (19.5%) 22.0 (11.1%)
    Non-smoker but previously smoked although not daily 163.0 (12.6%) 12.0 (5.9%) 19.0 (9.6%) 23.0 (10.5%) 32.0 (15.9%) 63.0 (23.2%) 14.0 (7.0%)
    Smoker but not daily 64.0 (4.9%) 6.0 (2.9%) 10.0 (5.1%) 7.0 (3.2%) 20.0 (10.0%) 12.0 (4.4%) 9.0 (4.5%)
    Unknown 3 0 0 1 2 0 0
Organic food (times/week) 0.5 (0.0, 3.0) 0.0 (0.0, 0.5) 0.5 (0.1, 3.0) 0.0 (0.0, 0.5) 1.0 (0.1, 3.0) 1.0 (0.5, 3.0) 0.0 (0.0, 1.0)
    Unknown 7 0 0 5 2 0 0
Pregnancy maternal active smoking 190.0 (15.1%) 25.0 (13.7%) 47.0 (23.7%) 55.0 (25.1%) 12.0 (6.0%) 9.0 (3.4%) 42.0 (21.2%)
    Unknown 40 22 0 2 4 11 1
Pregnancy maternal passive smoking 514.0 (40.3%) 55.0 (27.5%) 43.0 (21.8%) 126.0 (57.8%) 97.0 (48.7%) 14.0 (5.3%) 179.0 (90.4%)
    Unknown 21 4 1 3 4 8 1
Vegetables (times/week) 6.5 (4.0, 10.0) 6.0 (4.0, 10.0) 8.3 (4.4, 11.0) 6.0 (3.0, 8.5) 6.0 (3.5, 8.5) 8.5 (6.0, 14.0) 6.5 (4.0, 10.0)
    Unknown 6 1 0 2 2 1 0
a Median (IQR); n (%)
Table S7: Participant characteristics, by cohort and overall (HELIX subcohort; 2013-2016).

Concentrations of the glucocorticosteroids

Characteristic Overall, N = 1,004a BIB, N = 154a EDEN, N = 137a INMA, N = 205a KANC, N = 180a MOBA, N = 200a RHEA, N = 128a
Glucocorticosteroid
A 4.3 (2.4, 8.2) 4.8 (2.8, 9.0) 5.1 (2.6, 9.1) 3.0 (1.6, 5.6) 3.8 (2.0, 7.3) 4.3 (2.7, 8.4) 5.9 (3.5, 14.9)
    Unknown 1 0 0 1 0 0 0
E 22.9 (13.1, 38.5) 25.7 (14.5, 41.4) 28.6 (14.1, 42.0) 17.1 (10.3, 27.4) 21.4 (12.0, 33.7) 23.3 (14.1, 38.1) 28.9 (19.3, 59.4)
F 5.5 (3.2, 9.5) 6.3 (4.2, 10.4) 7.8 (4.2, 11.4) 4.6 (2.9, 7.1) 4.9 (2.7, 8.2) 5.2 (3.0, 9.1) 6.2 (3.4, 13.1)
    Unknown 2 0 0 0 1 1 0
Glucocorticosteroid metabolite
11OHAndros 234.2 (130.3, 390.5) 259.7 (151.9, 375.0) 413.0 (221.7, 617.0) 256.7 (142.9, 365.1) 163.3 (80.7, 298.5) 254.4 (151.5, 408.4) 165.4 (95.9, 304.2)
    Unknown 3 0 0 0 3 0 0
17-DO-cortolone 57.5 (29.1, 101.7) 56.1 (32.8, 100.6) 76.5 (46.0, 137.6) 61.3 (32.5, 102.1) 43.7 (15.1, 93.4) 56.4 (26.4, 92.0) 51.2 (28.5, 94.3)
    Unknown 2 0 0 0 1 0 1
20aDHE 16.6 (9.7, 27.5) 14.2 (7.0, 25.8) 25.8 (15.1, 37.8) 15.6 (10.2, 23.0) 14.8 (7.7, 25.6) 17.5 (11.7, 26.1) 14.8 (8.7, 27.6)
    Unknown 11 7 0 0 4 0 0
20aDHF 6.6 (3.3, 13.3) 7.2 (3.8, 14.0) 10.0 (5.7, 19.5) 5.5 (3.0, 9.4) 4.8 (2.2, 11.4) 7.4 (4.2, 14.0) 6.5 (2.9, 13.8)
    Unknown 7 4 0 0 3 0 0
20bDHE 9.5 (6.2, 14.3) 8.7 (4.8, 14.8) 13.2 (9.7, 17.3) 9.0 (6.6, 11.7) 8.9 (5.1, 13.7) 9.0 (5.9, 14.3) 8.7 (5.3, 15.2)
    Unknown 17 14 0 0 3 0 0
20bDHF 15.2 (9.1, 24.8) 16.5 (10.8, 26.5) 19.9 (12.0, 32.0) 13.0 (8.0, 18.1) 14.0 (8.5, 24.5) 14.2 (8.4, 23.5) 14.3 (7.9, 27.5)
5a20acortol 88.9 (52.1, 141.6) 109.8 (61.7, 177.3) 103.0 (58.0, 153.8) 83.0 (45.9, 118.7) 84.7 (46.9, 145.9) 88.6 (53.7, 138.2) 72.4 (47.2, 130.2)
    Unknown 9 9 0 0 0 0 0
5a20bcortol 122.4 (70.4, 185.0) 131.0 (66.3, 182.3) 148.8 (108.8, 226.1) 124.3 (68.9, 178.8) 115.2 (62.9, 189.2) 114.7 (67.8, 172.7) 105.3 (72.6, 175.0)
    Unknown 5 5 0 0 0 0 0
5aTHB 133.1 (76.1, 222.4) 159.8 (101.7, 241.3) 144.2 (87.9, 255.3) 115.7 (73.3, 171.7) 148.0 (82.6, 245.6) 106.1 (61.1, 184.9) 139.9 (74.6, 260.5)
5aTHE 73.9 (39.7, 124.0) 82.0 (52.1, 145.7) 83.9 (41.5, 132.7) 62.2 (32.3, 97.3) 71.3 (40.3, 121.7) 64.5 (36.4, 103.9) 107.9 (51.2, 183.2)
    Unknown 1 0 0 0 0 0 1
5aTHF 2,870.0 (1,663.7, 4,389.0) 3,394.6 (2,288.1, 5,308.1) 3,474.2 (1,856.1, 5,253.4) 2,756.9 (1,565.6, 3,758.3) 2,907.3 (1,656.1, 4,621.2) 2,283.3 (1,259.8, 3,454.6) 3,001.9 (1,652.3, 4,613.6)
5b20acortol 147.7 (83.5, 225.8) 177.4 (98.9, 302.3) 169.7 (91.1, 252.9) 141.9 (76.6, 187.6) 143.0 (80.2, 229.8) 143.7 (86.6, 204.2) 137.7 (79.6, 220.5)
    Unknown 11 11 0 0 0 0 0
5b20acortolone 641.9 (366.0, 983.1) 638.3 (385.0, 1,028.2) 903.7 (574.5, 1,296.1) 654.6 (398.7, 890.7) 518.0 (261.2, 870.2) 580.6 (318.0, 901.5) 629.3 (400.9, 962.4)
5b20bcortol 195.7 (120.1, 302.4) 242.7 (152.0, 356.8) 225.2 (142.1, 371.5) 199.9 (130.5, 289.3) 155.8 (88.0, 270.4) 186.3 (115.5, 269.4) 177.5 (113.7, 301.7)
    Unknown 3 3 0 0 0 0 0
5b20bcortolone 546.9 (336.3, 837.1) 561.3 (331.3, 889.9) 682.3 (452.0, 1,031.1) 534.1 (372.6, 792.7) 505.0 (272.3, 769.3) 496.1 (289.2, 761.3) 563.5 (328.4, 881.5)
5bDHF 1.4 (0.9, 2.0) 1.4 (0.9, 2.2) 1.8 (1.3, 2.6) 1.1 (0.6, 1.8) 1.5 (1.1, 1.9) 1.1 (0.6, 1.7) 1.5 (1.0, 2.1)
    Unknown 2 0 0 1 0 1 0
5bTHB 49.3 (28.0, 82.7) 53.3 (27.5, 98.3) 60.9 (34.9, 94.5) 50.0 (29.7, 73.1) 43.8 (27.5, 89.7) 40.0 (24.7, 65.7) 53.5 (28.4, 76.7)
    Unknown 1 0 0 0 1 0 0
5bTHE 3,138.3 (1,889.5, 4,694.0) 3,552.8 (2,335.3, 4,797.4) 3,649.6 (2,293.5, 5,317.1) 2,911.6 (1,615.2, 4,050.7) 2,754.6 (1,448.0, 3,989.3) 3,070.1 (1,785.5, 4,637.7) 3,541.6 (2,010.1, 5,901.3)
5bTHF 906.5 (548.0, 1,416.1) 1,116.2 (660.8, 1,644.8) 1,238.6 (743.1, 1,578.3) 882.9 (542.6, 1,199.8) 753.9 (389.4, 1,258.7) 859.7 (492.9, 1,261.3) 881.5 (565.0, 1,441.1)
    Unknown 2 2 0 0 0 0 0
6OHE 11.9 (6.5, 18.4) 13.2 (7.6, 20.6) 12.2 (6.1, 17.4) 9.2 (5.3, 14.1) 13.1 (7.1, 19.6) 11.2 (6.4, 18.1) 14.3 (8.7, 24.3)
6OHF 42.8 (22.5, 76.7) 51.9 (29.8, 93.9) 55.8 (29.8, 82.3) 32.3 (18.5, 53.3) 36.6 (19.7, 68.7) 46.0 (27.9, 82.9) 42.0 (21.1, 93.2)
Glucocorticosteroid precursor
S 0.4 (0.3, 0.8) 0.5 (0.3, 0.9) 0.4 (0.3, 0.7) 0.6 (0.4, 0.9) 0.3 (0.2, 0.5) 0.4 (0.3, 0.7) 0.4 (0.2, 0.8)
    Unknown 94 6 5 12 9 51 11
Glucocorticosteroid precursor metabolite
17HP 22.3 (15.1, 33.5) 17.0 (11.1, 27.6) 33.2 (23.5, 44.0) 20.3 (13.2, 32.2) 20.3 (10.8, 33.1) 23.0 (17.5, 31.2) 21.8 (15.7, 32.2)
    Unknown 1 0 0 0 0 0 1
5bDHS 0.3 (0.2, 0.4) 0.3 (0.2, 0.4) 0.3 (0.2, 0.5) 0.3 (0.2, 0.3) 0.2 (0.2, 0.3) 0.3 (0.2, 0.4) 0.3 (0.2, 0.5)
    Unknown 132 5 20 43 0 57 7
5bTHS 30.7 (18.5, 50.5) 35.7 (20.7, 59.2) 34.5 (19.8, 52.1) 27.7 (17.6, 43.0) 31.3 (18.6, 55.1) 26.2 (14.2, 40.8) 33.7 (20.0, 58.2)
    Unknown 2 0 0 1 0 1 0
PT 200.6 (112.8, 342.0) 149.1 (87.6, 246.3) 378.8 (230.8, 542.8) 253.4 (150.0, 404.4) 142.2 (82.4, 273.7) 176.4 (112.9, 283.3) 189.4 (104.9, 306.3)
Androgen
AED 0.2 (0.2, 0.3) 0.2 (0.2, 0.3) 0.3 (0.2, 0.5) 0.2 (0.1, 0.4) 0.2 (0.1, 0.3) 0.2 (0.1, 0.3) 0.2 (0.1, 1.1)
    Unknown 407 0 34 73 117 106 77
T 0.5 (0.3, 1.0) 0.7 (0.5, 1.0) 1.0 (0.5, 1.9) 0.6 (0.3, 1.0) 0.3 (0.2, 0.6) 0.4 (0.3, 0.7) 0.4 (0.3, 0.7)
    Unknown 75 0 5 3 29 24 14
Androgen metabolite
Andros 186.0 (78.1, 394.0) 148.4 (72.0, 267.9) 552.2 (308.7, 980.2) 295.4 (129.1, 513.8) 98.4 (39.6, 227.5) 134.7 (63.4, 293.1) 110.0 (61.6, 226.5)
    Unknown 1 0 0 0 1 0 0
Etio 110.9 (50.7, 237.8) 75.1 (32.6, 151.0) 369.7 (231.8, 561.0) 169.7 (84.0, 306.1) 74.8 (37.6, 122.6) 91.4 (45.8, 184.0) 76.2 (41.2, 147.0)
    Unknown 1 0 0 0 1 0 0
a Median (IQR)
Table S8: Participants glucocorticosteroids concentrations, by cohort and overall (HELIX subcohort; 2013-2016).

Tables for main results

Balancing weights: sample sizes

Exposure Unadjusted Adjusteda
Phenols
PRPA 1,297 1,297
ETPA 1,297 1,289
OXBE 1,297 1,277
BUPA 1,297 1,276
MEPA 1,297 1,266
TRCS 1,297 1,255
BPA 1,297 1,137
OP pesticide metabolites
DETP 1,297 1,222
DEP 1,297 1,222
DMTP 1,297 1,219
DMP 1,297 1,172
Phthalate metabolites
oxo-MiNP 1,297 1,199
oh-MiNP 1,297 1,171
MBzP 1,297 1,114
MEHP 1,297 1,090
MEP 1,297 1,054
MnBP 1,297 1,035
MEHHP 1,297 1,010
MEOHP 1,297 1,000
MECPP 1,297 980.4
MiBP 1,297 927.3
a Truncated weights.
Table S9: Effective sample size before and after balancing weights estimation (exposures: endocrine disrupting chemicals (EDCs); outcome: hit reaction time standard error (HRT-SE)) (HELIX subcohort; 2013-2016).
Exposure Unadjusted Adjusteda
Phenols
OXBE 976.0 960.1
PRPA 976.0 956.0
MEPA 976.0 953.7
BUPA 976.0 952.3
ETPA 976.0 951.7
TRCS 976.0 942.4
BPA 976.0 856.4
OP pesticide metabolites
DEP 976.0 922.1
DETP 976.0 922.1
DMTP 976.0 907.3
DMP 976.0 893.3
Phthalate metabolites
oh-MiNP 976.0 877.9
oxo-MiNP 976.0 873.6
MBzP 976.0 828.8
MEHP 976.0 827.3
MEP 976.0 796.3
MEHHP 976.0 784.8
MECPP 976.0 768.1
MEOHP 976.0 761.5
MnBP 976.0 745.7
MiBP 976.0 690.9
a Truncated weights.
Table S10: Effective sample size before and after balancing weights estimation (exposures: endocrine disrupting chemicals (EDCs); outcomes: glucocorticosteroids) (HELIX subcohort; 2013-2016).
Exposure Unadjusted Adjusteda
cortisone production 976.0 777.2
corticosterone production 976.0 757.5
cortisol production 976.0 751.5
a Truncated weights.
Table S11: Effective sample size before and after balancing weights estimation (exposures: glucocorticosteroids; outcome: hit reaction time standard error (HRT-SE)) (HELIX subcohort; 2013-2016).

Balancing weights: summary statistics

Characteristica Median (IQR) Range
N = 1,297a N = 1,297a
OP pesticide metabolites
DMP 0.99 (0.73, 1.25) 0.49, 1.50
DMTP 1.00 (0.81, 1.20) 0.59, 1.39
DEP 1.01 (0.81, 1.19) 0.59, 1.39
DETP 0.99 (0.81, 1.18) 0.61, 1.41
Phenols
MEPA 1.01 (0.90, 1.13) 0.74, 1.25
ETPA 1.01 (0.96, 1.07) 0.88, 1.14
PRPA

    2143289344 1,297 (100%) 1,297 (100%)
BPA 0.99 (0.70, 1.27) 0.39, 1.57
BUPA 1.01 (0.91, 1.11) 0.81, 1.22
OXBE 1.01 (0.92, 1.09) 0.79, 1.21
TRCS 1.01 (0.87, 1.13) 0.68, 1.28
Phthalate metabolites
MEP 0.93 (0.61, 1.27) 0.27, 1.77
MiBP 0.91 (0.46, 1.38) 0.05, 1.92
MnBP 0.98 (0.59, 1.33) 0.20, 1.74
MBzP 0.98 (0.66, 1.27) 0.35, 1.62
MEHP 0.98 (0.64, 1.28) 0.31, 1.68
MEHHP 0.96 (0.54, 1.35) 0.16, 1.76
MEOHP 0.96 (0.52, 1.35) 0.15, 1.78
MECPP 0.95 (0.50, 1.34) 0.14, 1.84
oh-MiNP 1.00 (0.74, 1.24) 0.47, 1.51
oxo-MiNP 1.01 (0.78, 1.20) 0.52, 1.43
a Truncated weights.
Table S12: Summary statistics of the estimated balancing weights (exposures: endocrine disrupting chemicals (EDCs); outcome: hit reaction time standard error (HRT-SE)) (HELIX subcohort; 2013-2016).
Characteristica Median (IQR) Range
N = 976a N = 976a
OP pesticide metabolites
DMP 0.99 (0.75, 1.23) 0.51, 1.46
DMTP 1.00 (0.78, 1.23) 0.56, 1.41
DEP 0.99 (0.81, 1.20) 0.64, 1.41
DETP 0.99 (0.82, 1.18) 0.62, 1.41
Phenols
MEPA 1.00 (0.90, 1.13) 0.75, 1.26
ETPA 1.02 (0.90, 1.14) 0.72, 1.24
PRPA 1.00 (0.92, 1.12) 0.76, 1.26
BPA 1.00 (0.70, 1.26) 0.40, 1.58
BUPA 1.01 (0.90, 1.13) 0.75, 1.27
OXBE 1.01 (0.92, 1.10) 0.78, 1.21
TRCS 1.01 (0.86, 1.14) 0.68, 1.29
Phthalate metabolites
MEP 0.92 (0.60, 1.27) 0.28, 1.74
MiBP 0.88 (0.44, 1.38) 0.09, 1.98
MnBP 0.97 (0.52, 1.35) 0.14, 1.84
MBzP 0.94 (0.68, 1.29) 0.35, 1.68
MEHP 0.98 (0.65, 1.29) 0.33, 1.64
MEHHP 0.98 (0.56, 1.35) 0.21, 1.69
MEOHP 0.98 (0.53, 1.35) 0.18, 1.77
MECPP 0.96 (0.55, 1.36) 0.19, 1.76
oh-MiNP 0.99 (0.73, 1.25) 0.45, 1.49
oxo-MiNP 1.01 (0.71, 1.25) 0.45, 1.52
a Truncated weights.
Table S13: Summary statistics of the estimated balancing weights (exposures: endocrine disrupting chemicals (EDCs); outcomes: glucocorticosteroids) (HELIX subcohort; 2013-2016).
Characteristica Median (IQR) Range
N = 976a N = 976a
cortisol production 1.00 (0.54, 1.39) 0.14, 1.80
cortisone production 1.00 (0.59, 1.39) 0.19, 1.73
corticosterone production 0.98 (0.56, 1.39) 0.15, 1.78
a Truncated weights.
Table S14: Summary statistics of the estimated balancing weights (exposures: glucocorticosteroids; outcome: hit reaction time standard error (HRT-SE)) (HELIX subcohort; 2013-2016).

Tables for other results

Balancing weights for effect modification: summary statistics

Characteristica Median (IQR) Range
females, N = 587a males, N = 710a females, N = 587a males, N = 710a
OP pesticide metabolites
DMP 0.99 (0.74, 1.25) 1.00 (0.74, 1.25) 0.53, 1.46 0.53, 1.46
DMTP 1.00 (0.79, 1.22) 1.01 (0.82, 1.20) 0.58, 1.38 0.58, 1.38
DEP 1.01 (0.82, 1.18) 1.02 (0.84, 1.17) 0.64, 1.36 0.64, 1.36
DETP 1.00 (0.77, 1.22) 1.01 (0.82, 1.20) 0.57, 1.39 0.57, 1.39
Phenols
MEPA 1.02 (0.89, 1.15) 1.02 (0.94, 1.11) 0.76, 1.23 0.76, 1.23
ETPA 1.02 (0.96, 1.08) 1.01 (0.97, 1.06) 0.91, 1.12 0.91, 1.12
PRPA 1.02 (0.92, 1.13) 1.02 (0.95, 1.10) 0.82, 1.21 0.82, 1.21
BPA 1.02 (0.73, 1.28) 1.02 (0.74, 1.25) 0.42, 1.50 0.42, 1.50
BUPA 1.02 (0.95, 1.10) 1.01 (0.81, 1.20) 0.67, 1.29 0.67, 1.29
OXBE 1.03 (0.92, 1.12) 1.02 (0.94, 1.09) 0.81, 1.19 0.81, 1.19
TRCS 1.03 (0.92, 1.13) 1.01 (0.89, 1.12) 0.73, 1.25 0.73, 1.25
Phthalate metabolites
MEP 0.96 (0.67, 1.26) 0.93 (0.62, 1.30) 0.31, 1.68 0.31, 1.68
MiBP 0.93 (0.51, 1.39) 0.96 (0.52, 1.40) 0.16, 1.85 0.16, 1.85
MnBP 1.00 (0.63, 1.33) 0.98 (0.59, 1.35) 0.28, 1.68 0.28, 1.68
MBzP 1.00 (0.71, 1.27) 0.99 (0.69, 1.27) 0.40, 1.57 0.40, 1.57
MEHP 1.02 (0.69, 1.27) 0.98 (0.62, 1.32) 0.33, 1.62 0.33, 1.62
MEHHP 1.01 (0.60, 1.29) 0.95 (0.56, 1.36) 0.26, 1.72 0.26, 1.72
MEOHP 1.00 (0.63, 1.29) 0.95 (0.53, 1.40) 0.23, 1.74 0.23, 1.74
MECPP 1.00 (0.59, 1.33) 0.95 (0.50, 1.37) 0.23, 1.76 0.23, 1.76
oh-MiNP 1.02 (0.78, 1.22) 1.00 (0.76, 1.23) 0.51, 1.46 0.51, 1.46
oxo-MiNP 1.02 (0.84, 1.17) 1.01 (0.76, 1.21) 0.58, 1.39 0.58, 1.39
a Truncated weights.
Table S15: Summary statistics of the estimated balancing weights for effect modification (exposures: endocrine disrupting chemicals (EDCs); outcome: hit reaction time standard error (HRT-SE); modifier: sex) (HELIX subcohort; 2013-2016).
Characteristica Median (IQR) Range
females, N = 434a males, N = 542a females, N = 434a males, N = 542a
OP pesticide metabolites
DMP 0.98 (0.77, 1.23) 1.01 (0.76, 1.21) 0.57, 1.45 0.57, 1.45
DMTP 1.03 (0.78, 1.22) 1.01 (0.79, 1.23) 0.56, 1.40 0.56, 1.40
DEP 1.01 (0.85, 1.16) 1.00 (0.84, 1.18) 0.67, 1.36 0.67, 1.36
DETP 1.00 (0.77, 1.22) 1.01 (0.86, 1.17) 0.57, 1.40 0.57, 1.40
Phenols
MEPA 1.01 (0.88, 1.16) 1.03 (0.94, 1.11) 0.73, 1.26 0.73, 1.26
ETPA 1.04 (0.92, 1.12) 1.02 (0.91, 1.12) 0.78, 1.22 0.78, 1.22
PRPA 1.03 (0.87, 1.16) 1.02 (0.95, 1.10) 0.74, 1.24 0.74, 1.24
BPA 1.00 (0.71, 1.28) 1.01 (0.75, 1.24) 0.44, 1.52 0.44, 1.52
BUPA 1.02 (0.95, 1.11) 1.01 (0.80, 1.20) 0.64, 1.30 0.64, 1.30
OXBE 1.03 (0.86, 1.16) 1.02 (0.95, 1.09) 0.76, 1.22 0.76, 1.22
TRCS 1.03 (0.92, 1.13) 1.01 (0.88, 1.14) 0.73, 1.25 0.73, 1.25
Phthalate metabolites
MEP 0.99 (0.70, 1.24) 0.95 (0.55, 1.30) 0.31, 1.68 0.31, 1.68
MiBP 0.92 (0.46, 1.40) 0.92 (0.54, 1.39) 0.15, 1.84 0.15, 1.84
MnBP 0.97 (0.51, 1.40) 0.98 (0.57, 1.32) 0.21, 1.78 0.21, 1.78
MBzP 0.99 (0.70, 1.26) 0.98 (0.66, 1.31) 0.38, 1.58 0.38, 1.58
MEHP 1.01 (0.72, 1.29) 0.98 (0.61, 1.34) 0.36, 1.58 0.36, 1.58
MEHHP 1.02 (0.65, 1.31) 1.00 (0.59, 1.35) 0.30, 1.63 0.30, 1.63
MEOHP 1.01 (0.62, 1.32) 1.01 (0.51, 1.41) 0.24, 1.68 0.24, 1.68
MECPP 0.98 (0.62, 1.32) 0.98 (0.54, 1.40) 0.29, 1.67 0.29, 1.67
oh-MiNP 1.00 (0.73, 1.26) 1.00 (0.78, 1.24) 0.49, 1.44 0.49, 1.44
oxo-MiNP 1.03 (0.74, 1.27) 1.02 (0.76, 1.24) 0.47, 1.45 0.47, 1.45
a Truncated weights.
Table S16: Summary statistics of the estimated balancing weights for effect modification (exposures: endocrine disrupting chemicals (EDCs); outcomes: glucocorticosteroids; modifier: sex) (HELIX subcohort; 2013-2016).
Characteristica Median (IQR) Range
females, N = 434a males, N = 542a females, N = 434a males, N = 542a
cortisol production 0.97 (0.57, 1.41) 1.01 (0.59, 1.35) 0.24, 1.71 0.24, 1.71
cortisone production 1.00 (0.61, 1.40) 1.00 (0.59, 1.38) 0.27, 1.69 0.27, 1.69
corticosterone production 1.00 (0.60, 1.39) 1.03 (0.56, 1.37) 0.23, 1.71 0.23, 1.71
a Truncated weights.
Table S17: Summary statistics of the estimated balancing weights for effect modification (exposures: glucocorticosteroids; outcome: hit reaction time standard error (HRT-SE); modifier: sex) (HELIX subcohort; 2013-2016).

Supplementary figures

Figures for descriptive data

Study populations

flowchart TB
  helixsc["HELIX subcohort\n(N = 1,301)"] --> edcs["HELIX data EDCs\n(N = 1,297)"]
  helixsc --> corts["HELIX data glucocorticosteroids\n(N = 1,004)"]
  edcs --> inter["HELIX data EDCs and glucocorticosteroids\n(N = 976)"]
  corts --> inter
  edcs -.-> rq1(["HRT-SE ~ EDCs"])
  inter -.-> rq2(["metabolites ~ EDCs"])
  inter -.-> rq3(["HRT-SE ~ metabolites"])

Figure S1: Flowchart describing the sample size for each research question.

Description of endocrine disruptors

Figure S2: Measurement classification of endocrine disrupting chemicals (EDCs), by cohort (HELIX subcohort; 2013-2016). Coding: 1, quantifiable; 2, <LOD; 3, interference or out of range; 4. not analysed.

Description of glucocorticosteroids

Figure S3: Measurement classification of the glucocorticosteroids, by cohort (HELIX subcohort; 2013-2016). Coding: 1, quantifiable; 2, <LOQ; 3, interference or out of range; 4, not detected.

Figures for other results

Marginal contrasts for effect modification

Figure S4: Marginal contrasts on the logarithmic scale for effect modification by sex of a increase from the 10th to the 90th percentile of the endocrine disrupting chemicals (EDCs) on hit reaction time standard error (HRT-SE) expressed in ms (HELIX subcohort; 2013-2016). Circles and triangles indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56).

Figure S5: Marginal contrasts on the logarithmic scale for effect modification by sex of a increase from the 10th to the 90th percentile of the endocrine disrupting chemicals (EDCs) on the glucocorticosteroids expressed in ng/ml (HELIX subcohort; 2013-2016). Circles and triangles indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56).

Figure S6: Marginal contrasts on the logarithmic scale for effect modification by sex of a increase from the 10th to the 90th percentile of the glucocorticosteroids on hit reaction time standard error (HRT-SE) expressed in ms (HELIX subcohort; 2013-2016). Circles and triangles indicate effect estimates. Solid lines indicate the \(95\%\) CI. The size of the circles represents the S value of the effect estimate (56). Abbreviations: cortisone production (cortisone prod.); cortisol production (cortisol prod.); corticost. prod. (corticosterone production).